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FOURIER ANALYSIS METHODS FOR THE COMPRESSIBLE 
NAVIER-STOKES EQUATIONS 

RAPHAEL DANCHIN 


Abstract. In the last three decades, Fourier analysis methods have known a growing 
importance in the study of linear and nonlinear PDE’s. In particular, techniques based 
on Littlewood-Paley decomposition and paradifferential calculus have proved to be very 
efficient for investigating evolutionary fluid mechanics equations in the whole space or in 
the torus. We here give an overview of results that we can get by Fourier analysis and 
paradifferential calculus, for the compressible Navier-Stokes equations. We focus on the 
Initial Value Problem in the case where the fluid domain is (or the torus with 
d > 2, and also establish some asymptotic properties of global small solutions. The time 
decay estimates in the critical regularity framework that are stated at the end of the 
survey are new, to the best of our knowledge. 


1. Introduction 

In the Eulerian description, a general compressible fluid evolving in some open set 11 
of is characterized at every material point x in H and time t € R by its velocity field 
u = u(t, x) E R*^, density g = g{t, x) E R+, pressure p = p{t, x) E R, internal energy by unit 
mass e = e{t, x) E R, entropy by unit mass s = s(t, x) and absolute temperature T = T(t, x). 
In the absence of external forces, those quantities are governed by: 

• The mass balance: 

dtg + div (gu) = 0. 

• The momentum balanc43: 

dt{gu) + div {gu ®u) = divr — Vp, 

where r stands for the viscous stress tensor. 

• The energy balance: 

dt(^Q{e + + div (^g{e + = div(T • u + pu) - divq, 

where q is the heat flux vector. 

• The entropy inequality: 

(1) dt{gs) + diY (gsu) >-div [§^)- 

In what follows, we concentrate on so-called Newtonian fluid]^. Hence (see e.g. m) T is 
given by: 

r'^^Adivttld +2pD{u), 

^With the convention that (div (o (8) 6))^ di(a* b^). 

^That is to say: the viscous stress tensor -r is a linear function of D^u, invariant under rigid transforms, 
there is no internal mass couples (and thus the angular momentum is conserved), and the fluid is isotropic 
(viz. the physical quantities depend only on {t,x)). 
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where the real numbers A and ^ are the viscosity coefficients and D{u) \ {Du + is 

the deformation tensor. 

If we assume in addition that the Fourier law q = —k\/T is satished then we get the 
following system of equations: 

{ dtQ + div(^u) = 0, 

dt{Qu) + div(^u (g) tt) — 2div (/rL)(n)) — V(Adivn) + Vp = 0, 
dtlge) + div(^eu) +pdivu — div(A:VT) = 2pD{u) : D{u) + A(divu)^. 

We shall further postulate that the entropy s is interrelated with p, T and e through the 
so-called Gibbs relation 

Tds = de+pd(^-^, 
and thus we get the following evolution equation for s: 

(3) T[dt{gs) + div(ijsu)) = r • D{u) — divg. 


For the entropy inequality to be satisfied, a necessary and sufficient condition is thus 

ivrP 

T : D{u) k ' > 0, 


which yields the following constraints on A, p and k: 


k > 0, p >0 and 2p + dX> 0. 


In order to close System ([2]) which is composed of d+2 equations for d+A unknowns (namely 
Q, e, p, T and ,■ ■ ■ , u'^), we need another two state equations interrelating p, g, e, s and T. 
In this survey, for simplicity we shall focus on barotropic gases that is p depends only on 
the density and A and p are independent of T. Therefore the system constituted by the 
first two equations in ([2]), the so-called barotropic compressible Navier-Stokes system : 


f 9i£»-hdiv(£»u) = 0, 

( dt{gu) + div Iqu 0 u) — div (2pD{u) + Adivuld) -|- Vp = 0 


where p 


def 


P{g) for some given smooth function P, is closed. 


Our main goal is to solve the Initial Value Problem (or Cauchy Problem) for ([4]) sup¬ 
plemented with initial data {go,uo) at time t = 0 in the case where the fluid domain fl 
is the whole space or the torus. We will concentrate on the local well-posedness issue for 
large data with no vacuum, on the global well-posedness issue for small perturbations of 
a constant stable equilibrium, and will give exhibit some of the qualitative properties of 
the constructed solutions. As regards global results, the concept of critical regularity is 
fundamental. Indeed, experience shows that whenever the PDF system under considera¬ 
tion possesses some scaling invariance with respect to space and time dilations (which is 
in general the case when it comes from mathematical physics) then appropriate so-called 
critical norms or quantities essentially control the (possible) finite time blow-up and the 
asymptotic properties of the solutions. 

If modifying the pressure law accordingly then the barotropic system (j3]) we are here 
considering has the following scaling invariance: 

(5) g{t,x) g{i‘^t,ix), u{t,x) iu{i^t,ix), i > 0. 
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More precisely, {g, u) is a solution to (jH) if and only if so does ui), with pressure function 
£^P. This means that we expect optimal solution spaces (and norms) for (j3|) to have the 
scaling invariance pointed out above. 

The rest of these notes unfolds as follows. In the next section, we present the basic 
tools and estimates that will be needed to study System ([261) . Then we concentrate on the 
local well-posedness issue for ([26|) in critical spaces. Section H] is dedicated to solving ([1|) 
globally for small data. The last section is devoted to asymptotic results for the system. 
We concentrate on the low Mach number limit and on the decay rates of global solutions 
in the critical regularity framework. 

2. The Fourier analysis toolbox 

We here shortly introduce the Fourier analysis tools needed in this survey, then state 
estimates for the heat and transport equations that will play a fundamental role. For the 
sake of conciseness, some proofs are just sketched or omitted. Unless otherwise specified, 
the reader will find details in [T], Chap. 2 or 3. 

2.1. The Littlewood-Paley decomposition. The Littlewood-Paley decomposition is a 
dyadic localization procedure in the frequency space for tempered distributions over 
One of the main motivations for using it when dealing with PDFs is that the derivatives 
act almost as dilations on distributions with Fourier transform supported in a ball or an 
annulus, as regards norms. This is exactly what is stated in the following proposition: 

Proposition 2.1 (Bernstein inequalities). Let 0 < r < R. 

• There exists a constant C so that, for any k G N, couple {p,q) in [l,oo]^ with 
q > p > 1 and function u of with u supported in the ball B{0,XR) o/ for 
some X > 0, we have D^u G L*? and 

• There exists a constant C so that for any k G N, p G [l,oo] and function u of 
with Supp u C G / rA < 1^1 < RX} for some X > 0, we have 

A^nlliP < C^+^\\D’^u\\lp. 

As general solutions to nonlinear PDF’s need not be spectrally localized in annuli, we 
want a device for splitting any function or distribution into a sum of spectrally localized 
functions. To this end, fix some smooth radial non increasing function y with Suppy C 
B{0, |) and x = 1 on B{0, |), then set = x(C/2) — x(0- We thus have 

= 1 in 
jez 

The homogeneous dyadic blocks Aj are defined by 

Aju '^=ip{2-W)u T-'^iip{2-^-)Pu) = 2^^h{2^-)iru with h '^= p-^Lp. 

We also introduce the low frequency cut-off operator Sj: 

Sju x(2"^D)u = 2^^h{2^-) A u with h p-\. 

Let us emphasize that operators Aj and Sj are continuous on L^, with norm independent 
of j, a property that would fail if taking a rough cut-off function x (unless p = 2 of course). 


4 


RAPHAEL DANCHIN 


The price to pay for smooth cut-off is that Aj is not an orthogonal projector. However 
the following important quasi-orthogonality property is fulfilled: 

(6) AjAk = 0 if \j - k\ > 1. 

The homogeneous Littlewood-Paley decomposition for u reads 

(7) '^ = '^ 

3 

This equality holds modulo polynomials only. In order to have equality in the distributional 
sense, one may consider the set 5^ of tempered distributions u such that 

lim ||5ju||Lcx) = 0. 
j^-oo 

As distributions of 5^ tend to 0 at inhnity, one can easily conclude that ([7]) holds true in 
S' whenever u is in 5^. 


2.2. Besov spaces. It is obvious that for all s G M, we have 

(8) < C||u||^,, 

jGZ 

and it is also not very difficult to prove that for s G (0,1), 

(9) C'“^||u||(^o,s < sup2'^^||Ajtt||Loo < C'||u||(j,o.s, s G (0,1). 

jez 

In ([8|) and Q, we observe that three parameters come into play: the regularity parameter 
s, the Lebesgue exponent that is used for bounding AjU and the type of summation that 
is done over Z. This motivates the following dehnition: 


Definition 2.1. For s G M and 1 < p,r < oo, we set 
def 


\u\\bs — 

p,r 


def 


y^2^'^'^||A,'n||rp ) if r < oo and ||u||rjs =''sup2 -^^||Aj-uIHp. 


ITe then define the homogeneous Besov space to he the subset of distributions u G S'f^ 
such that ||u||ns < oo. 

^p,r 

We shall often use the following classical properties: 

• Scaling invariance'. For any s G K and (p, r) G [l,-|-oo]^ there exists a constant C 
such that for all A > 0 and u G Bp^,, we have 


(10) 


C p ||u||< ||^t(A- 


|o„ < CA^ ^IkllR^ • 

^p,v ^p,r 


Completeness'. Bp^. is a Banach space whenever s<d/p or s<d/p and r = 1. 
Fatou property: if {un)neN is a bounded sequence of functions of ,, that converges 
in S' to some u G S'f^ then u G Bf ,. and ||u||^s < C lim inf ||u^ 

Duality: If u is in then we have 

ll^ill^s <Csup{u,(j)) 




where the supremum is taken over those (/> in 5 n B such that ||(/>||^-s < 1. 
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Interpolation: The following inequalities are satisfied foiH all 1 < p,Fi,r 2 ,r < oo, 
Si / S 2 and 0 E (0,1): 


I 11 11 11 ^ _ 0 


\U\ 




d»2 ■ 

JPp,r2 


Action of Fourier multipliers: If T is a smooth homogeneous of degree m function 
on \ {0} then 


( 11 ) 


F{D) : k 


S 

p,r 


Tbs—m 


In particular, the gradient operator maps ,, in B, 


s-l 


Proposition 2.2 (Embedding for Besov spaces on M'^). (1) For any p E [l,oo] we have 

the continuous embedding i ^ IP ^ -6p,oo' 

. s—d(— -i-) 

(2) //s E M, 1 < Pi < p 2 < oo and 1 < ri < r 2 < oo, then ^ ^P 2 ,r 2 ^^ ■ 

(3) For s' < s and any 1 <p,ri,r 2 < oo, the embedding ofBp j.^ in Bp j.^ is locally com- 

. s'—d{— -5-) 

pact, i.e. for any (p £ S, the map pu is compact from to Bp 2 ,r 2 • 

. ^ 

(4) The space Bp^ is continuously embedded in the set of hounded continuous functions 
(going to 0 at infinity if, additionally, p < oo). 

2.3. Paraproduct and nonlinear estimates. Formally, the product of two tempered 
distributions u and v may be decomposed into 

(12) uv = TuV + R{u,v)+T„u 

with 

TuV ^Sj-iu AjV and B(u, AjuAj/v. 

i i 

The above operator T is called “paraproduct” whereas R is called “remainder”. The 
decomposition (jl2p has been first introduced by J.-M. Bony in [3]. We observe that in 
Fourier variables the sum in T^v is locally finite, hence is always defined. We shall see 
however that it cannot be smoother than what is given by high frequencies, namely v. As 
for the remainder, it may be not defined, but if it is then the regularity exponents add up. 
All that is detailed below: 

Proposition 2.3. Let (s, r) E M x [1, oo] and 1 < p,pi,P 2 < cc with 1/p = 1/pi + l/p 2 - 
• We have: 


\T,v 


Irs ^ ll^ll Rs and llruPlI/js+t ^ II’^IIr* I 

^p,r °P2,^ ^P:^ ^Pl,oo 


, if t <0. 


If Si + S 2 > 0 and l/r = l/ri + l/r 2 <l then 


||B(u,n)||psi+=2 < II'wIIru lbllR=2 . 


// Si + S 2 = 0 and l/ri + l/r 2 >l then 


||B(n,u)||, 


W,v}\\^o < ||n| 

P5OO 


^Pl,'ri -^P2’^2 


Putting together decomposition (I12h and the above results, one may get the following 
product estimate that depends only linearly on the highest norm of u and v: 


^With the convention that A < B means that A < CB for some ‘harmless’ positive constant C. 
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Corollary 2.1. Let u and v be in n Sp ,, for some s > 0 and {p,r) € [l,oo]^. Then 
there exists a constant C depending only on d, p and s and such that 


uv 


<C{\\ 




+ IhllLoo \\u\\ 


)■ 


Remark 2.1. Because is embedded in L°°, we deduce that whenever p < +oo, the 

d d 

product of two functions in B^^ is also in B^^ and that for some eonstant C = C{p,d): 

Haft’ll d < C'll'^^ll d Ilf II d . 

BP, BP, BP, 

P,1 p,l p,l 

Let us finally state a composition result. 


Proposition 2.4. Let F : M — >• M 5e smooth with F(0) = 0. For all 1 < p,r < oo and 
s > 0 we have F{u) G ,, n L°° for u G n L°°, and 

(13) l|i"(«)llB. <C\\u\\j,, 

p,r 

with C depending only on ||u||i;,oo, F' (and higher derivatives), s, p and d. 


2.4. Endpoint maximal regnlarity for the linear heat equation. This paragraph is 
dedicated to maximal regularity issues for the basic heat equation 

(14) dtu - An = /, n|t=o = uq. 

In the case no = 0, we say that the functional space X endowed with norm || • ||x has the 
maximal regularity property if 

(15) \\dtu,Dlu\\x<C\\f\\x. 

Fourier-Plancherel theorem implies that (|15p holds true for X = x M'^). From more 

complicated tools based on singular integrals and heat kernel estimates, one may gather 
that (jl5j) is true for X = L^(R+; LP(]R'^)) whenever 1 < p,r < +oo. However, the endpoint 
cases where one of the exponents p or r is 1 or +oo are false. 

One of the keys to the approach presented in these notes is that (fT5]l is true for X = 
L^(M+; Bp ]^(M'^)), a consequence of the following lemma: 

Lemma 2.1. There exist two positive constants c and C such that for any j £ Z, p ^ [1, oo] 
and A G R+, we have for all u £ S' with AjU in L^ : 

\\e^^Aju\\LP < llAjnlliP. 

Proof. A suitable change of variable reduces the proof to the case j = 0. Then consider a 
function (j) in \ {0}) with value 1 on a neighborhood of Suppp. We have 

= gx-kAou with g'A(a:) (27r)“'^ ^ 

Integrating by parts yields 

gx{x) = (1 + |x|2)-'^^^F("l«)(Id - 

Combining Leibniz and Faa-di-Bruno’s formulae to bound the integrant, we get 

|5a(x))| <C(1 + N2)-V‘=°\ 


FOURIER ANALYSIS METHODS FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS 


7 


and thus 

(16) IIsaIIli < 

Now the desired inequality (with j = 0) just follows from —)• LP. 


□ 


Theorem 2.1. Let u satisfy m- Then for any p G [l,oo] and s € M the following 
inequality holds true for all t > 0 : 


(17) 


k(i)llj 


||V^m||^s ^ dr < C^lluoll^a ^ + 


3“ , 
p,i 


dr 


Proof. If u satisfies (1141) then we have for any j G Z, 

Aju{t) = e^^AjUo + f Aj f (t) dr. 

Jo 

Taking advantage of Lemma 12.11 we thus have 

(18) ||A,n(t)||rP <e-'=°2^^'*||A,no||LP+ T e-'^°2'^'(*-")||A,/(r)||iP dr. 

Jo 

Multiplying by 2^^ and summing up over j yields 


J 

whence 




t)||lp dr 


\u\ 


~ -r \\J \\Li{0,t-,B‘ p- 


~ ll'U'O Ilijs + 
p.i 

Note that integrating (fTSll with respect to time also yields 

2^1||Ajtt|Ai(o,t;LP) ^ (l - e (^WAjUoWiP + \\Ajf\\L'L{o,t-,Lp)J- 

Therefore, multiplying by 2^^, using Bernstein inequality and summing up over j yields 
(19) \\^MLH0,t;B;^p ^ E(l - e-^°^"'^)2^^[\\Ajuo\\LP + ||A,/|Ai(o,t;rP)), 


which is even slightly better than what we wanted to prove 


□ 


Remark 2.2. Starting from (1181) and using general eonvolution inequalities in 1R+ gives 
a whole family of estimates for the heat equation. However, as time integration has been 
performed before summation over j, the norms that naturally appear are 


def 




2JC 




where 






def I 


Il“((0, t);L*> («<*))• 


With this notation, dUD implies that 


IkL ^Ikollijs +II/IL„ .s-2+^ for l<p2<Pi<oo. 

The relevancy of the above norms in the maximal regularity estimates has been first noticed 
(in a particular case) in the pioneering work by J.-Y. Chemin and N. Lerner [6], then 
extended to general Besov spaces in [5] . They will play a fundamental role in the proof of 
decay estimates, at the end of the paper. 

Let us point out that results in the spirit of Propositions\2fM and \2.4\ may be easily proved 
for LKBp j.) spaces, the general rule being just that the time exponents behave according to 
Holder inequality. 
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2.5. The linear transport equation. Here we give estimates in Besov spaces for the 
following transport equation: 


( 20 ) 


in 

in 


dta + V ■ Vxa + Ao = / 
a\t=o = 

where the initial data oq = ao{x), the source term / = f{t,x), the damping coefficient 
A > 0 and the time dependent transport field v = v{t,x) are given. 

Assuming that uq £ X and / E A), the relevant assumptions on v for (f20]l to 

be uniquely solvable depend on the nature of the Banach space X. Broadly speaking, in 
the classical theory based on Cauchy-Lipschitz theorem, v has to be at least integrable in 
time with values in the set of Lipschitz functions, so that it has a flow i/^. This allows to 
get the following explicit solution for ([20|) : 

(21) a(t,x) = e"^*ao(V’t“^(a:)) + / dr. 

JQ 

Theorem 2.2. Let 1 < p < < oo, 1 < r < oo and s E M satisfy 

. f d d\ , d 

— mm —,— <s<lH - 

\pi p'J Pi 

Then any smooth enough solution to (l20l) fulfills 


lallrcxj/6., ^ + A||a||^i^^^< e ^ 




“oll^s + 




with 


V{t) if 


I|Vi>(t)|| i dr. 


10 


rPI n^o 

-^Pl ,oo' 


In the ease s = 1 + — and r = 1, the above inequality is true with V'{t) = ||Vn(t)|| _d_ . 

Pi ^pi 

pi.i 


Proof. Applying Aj to ([2011 gives 

(22) dfAjO + V ■ VAjO + AAjO = Ajf + Rj with Rj [n • V, Aj]a. 

Therefore, from classical estimates for the transport equation, we get 


def 


( 23 ) ||Aja(t)||i;,p + A||Aja|Ai(^p) < ||Ajao||i,p 

+ ^*(||A,/|Ap + WRjWlp + ||A,a|Ap) dr. 

We claim that the remainder term Rj satisfies 

( 24 ) \\Rj{t)\\LP < Ccj{t) 2 ~f‘'\\Xv{t)\\ jl ||a(i)|lB» \\{cj{t))\\e-= I- 

4A.ocnL- 

Indeed, from Bony’s decomposition, we infer that (with the summation convention over 
repeated indices): 

( 25 ) Rj = [T,, , Aj]dka + - AjPg^y + i^(n^ A,a) - AjR{v\dka). 

The first term is the only one where having a commutator improves the estimate. Indeed, 
owing to the properties of spectral localization of the Littlewood-Paley decomposition, we 
have 

[T^k,Aj]dka= ^ [Sj>_iv^,Aj]dkAj,a. 
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Now, remark that 

[Sy_iv'^,Aj]dkAya{x) = 2^'^ [ h{2^(x - y)){Sy_iv'‘{x) - Sj>_iv'‘{y))dkAj>a{y) dy. 
Hence using the mean value formula and Bernstein inequalities yields 

||[r,.,A,]4a||LP <2-^'||5,._iVu||ioc ^ \\dkAj,a\\LP<\\Vv\\Loo ^ \\Aj,a\\LP. 

Bounding the third and last term in (|25p follows from Proposition 12.31 Regarding the 
second term, we may write 




Sji-idkAja A 




hence using Bernstein inequality, 


~ 2^ ^ ||Aja||LP ||vAj/u||l-, 

f>j-s 

and convolution inequality for series thus ensures (j24p for that term. 

Finally, we have 

AjR{v^,dka) = dkAjia (Aj/_i + Af + Aj/+i)u^, 

whence, by virtue of Bernstein inequality, 

\\R{v^,dka)\\LP < \\Aj>a\\Lp\\'^v\\L-o, 

and that term is thus also bounded by the r.h.s. of (1241) . 

Let us resume to (j22p . Using (1231) and (1240 . multiplying by 2^^ then summing up over j 
and using the notations of Remark 12.21 yields 

Then applying GronwalPs lemma gives the desired inequality for a. 


I n 11 jhs d'T . 

^p,r 


□ 


3. The local existence in critical spaces 


This section is dedicated to solving ([3P locally in time, in critical spaces. For simplicity, 

we focus on the case where the density goes to 1 at oo. Setting a — 1 and looking for 
reasonably smooth solutions with positive density. System (j4j) is equivalent to 

{ dta + u ■ Va = —(1 + a)divri, 

dtu + u ■ Vu- — -h VG(a) =- div {2'jl(a)D(u) + A(a)divu Id), 

1 + a 1 + a 

wher^ G'{a) ^ + (A + ^)Vdiv with A A(l) and y Ai(l), 

'jl{z) y,{l + z) — fi{l) and A(z)A(1 + z) — A(l). 


'*In our approadi, tlie exact value of functions G, A and y will not matter. We shall only need enough 
smoothness, and vanishing at 0 for A and Jl. 
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The scaling invariance properties for {a,u) are those pointed out for (jH). Critical norms 
for the initial data are thus invariant for all £ > 0 by 

ao{x) ao{£x) and uo{x) iuo{ix). 

In all that follows, we shall only consider homogeneous Besov spaces having the above 
scaling invariance and last index 1. There are good reasons for that choice, which will be 
explained throughout. Remembering (jlOl) . we thus take 

(27) ao € Bp I and uq ^ • 

In order to guess what is the relevant solution space, we just use the fact that a is governed 
by a transport equation and that u may be seen as the solution to the following Lame 
equation: 


(28) dtu - Au = f, u\t=o = UQ. 

In the whole space case, the solutions to (I28p also satisfy the estimates of Theorem 12. 1 1 and 
Remark 12.21 whenever the following ellipticity condition is fulfilled: 

(29) iJ.>0 and i/ = A + 2/r > 0. 

Indeed, if we denote by V and Q the orthogonal projectors over divergence-free and poten¬ 
tial vector fields, then we have 

dtVu — fiAVu = V f and dtQu — uAQu = Qf. 

In particular, applying Theorem 12.11 yields for all t >0: 


\\Vu{t)\\ 

WQums 


pA 


ill ||V^Pu||^. dr <c(\\Vuo\\^s 
Jo P’l V 




\\Vf\\ 


B. 


dr 


pA 


pA 


+ u \\V^Qu\\^s dr <(7* II QuoII os -b 


P-1 


P-1 




dr 


As V and Q are continuous on (being 0 order multipliers) we conclude that 


(30) ||'w(i)|lB'* -b min 


P-1 


Jo 


V^u||os dr < Cl ||uo|| 


P-1 


P-1 


/‘ 

Jo 


B. 


dr 


p,i 


• ^-1 

So, in short, starting from uo G Bp^ , we expect u in (|26D to be in 

Ep{T) {u G C([0,r];R|”\ dtU,V^u G L^O, T; r|"')}- 
In particular Vn has exactly the regularity needed in Theorem 12.21 to ensure the conserva- 

d 

tion of the initial regularity for a, and we thus expect a G C([0, T]; Bp^). 

The rest of this section is devoted to the proof of the following statement: 

Theorem 3.1. Let the viscosity coefficients A and // depend smoothly on g, and satisfy 
([2^ . Assume that (ao,uo) fulfills (127)1 for some 1 <p < 2d, and that d > 2. If in addition 

1-bao is bounded away from 0 then (l2^ has a unique local-in-time solutioi^ and (a,u) with 
. 

a in C{[0,T]; Bp anduinEp{T). 

^Owing to Remark 12.21 and to Theorem [U the constructed solution will have the additional property 
that u € ) and that a G 
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We propose two different proofs for Theorem 13.11 The first one uses an iterative scheme 
for approximating (I26p which consists in solving a linear transport equation and the Lame 
equation with appropriate right-hand sides. Taking advantage of Theorem 12.21 and of (f30|) . 
it is easy to bound the sequence in the expected solution space on some fixed time interval 
[0, T] with T > 0. However, because the whole system is not fully parabolic, the strong 
convergence of the sequence is shown for a weaker norm corresponding to ‘a loss of one 
derivative’. For that reason, that approach works onl}0 if 1 < p < d and d > 3. The 
same restriction occurs as regards the uniqueness issue, although the limit case p = d, or 
d = 2 and p < 2 is tractable by taking advantage of a logarithmic interpolation inequality 
combined with Osgood lemma (see the end of this section). 

The second proof consists in rewriting (I26p in Lagrangian coordinates. Then the density 
becomes essentially time independent, and one just has to concentrate on the velocity 
equation that is of parabolic type for small enough time, and can thus be solved by the 
contracting mapping argument. 

3.1. The classical proof in Eulerian coordinates. We here present the direct approach 
for solving ([26|) . Our proof covers only the case d>3 and 1 < p < d as regards existence, 
and 1 < p < d with d > 2 for uniqueness (variations on the method would allow to get 
existence for the full range 1 < p < 2d with d > 2, though). To simplify the presentation, 
we assume that A and p are density independent so that (I26|] rewrites 


dta + u ■ Va = —(1 + a)divu, 

dtu — Au = —u ■ Vn — I{a)Au — V(G(a)) 



Furthermore, we suppose that for a small enough constant c = c{p, d, G), 


(31) 


d < c. 


Step 1: An iterative scheme. We set Og = Snao and Uq = SnUo, and define the first term 
(a'^,n'^) of the sequence of approximate solutions to be 



where (e*'^)t>o stands for the semi-group of operators associated to (1281) . 


Next, once has been constructed, we define and to be the solutions 

to the following linear transport and Lame equations: 



supplemented with initial data Og^^ and Ug^^ 


®Let us emphasize however that one may modify the iterative scheme for constructing solutions, then 
resort to compactness arguments to get existence for the full range 1 < p < 2d and d > 2. 
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Step 2: Uniform estimates in the case 1 < p < 2d and d > 2. As the data are smooth, 
it is not difficult to check (by induction) that and are smooth and globally defined. 

We claim that there exists some T > 0 such that {a"')n£n is bounded in C{[0,T]; and 

’ _ d 

{u^)n£N is bounded in the space Ep{T). Indeed, Theorem 12.21 and the fact that is 
stable by product imply that for some C > 1, 

rT 

\a^^\T)\\ d < ||a"+i| 


RP, 


\ d +C [ (l + ||a"'|| d )||divu"'|| d dt 


+C 


||Vn^ 


d lla”’'"^ll d dt. 

P 
P^ 


BP, 




Let U^{T) ||Vn”'|| d . Applying Gronwall’s lemma and using the definition of 

we thus get 


,CU^{T) 


- 1 


(33) d < d + (1 + ||a^dl i 

Therefore, assuming that oq fulfills (j3ip with some small enough c, that 

(34) lla^ll d <4Cc 

and that 

(35) (7C/”(r) <log(l + c), 

we conclude from ([33|) that also satisfies (fM)) for the same T. At this point, let us 

d 

observe that as B^, is continuously embedded in L°°, one may take c so small as 

(36) ||o"'|| d <4Cc implies ||a"'||LcxifroTlxK'i'i ^ 1/2- 

Let us now prove estimates for the velocity. From (|3UD . we get for some constant C 
depending only on A and /i, 

rT 


W 


n+l I 


EJT)<c(\\uo\\ d_,+ [ \\u'^-W + I{a^)Au^ + ViG{c 

V Jo 


I d_. dt 

BP, 

P,1 


The terms in the r.h.s. may be bounded by means of Propositions 12. 3l and 12. dl lremembering 
(l36l) l if d > 2 and 1 < p < 2d. We get for some C = C"(p, d, G): 


n+l 


Irp(T) < Ikoll . j.-i + (II 

A B , 


a'l\ d +||u"|| d_i)||Vu"|| d 


L¥iBfi) 


B^iBU 1' 


+r||a" 




Using ()34p and the definition of this implies that 


u 


n+l I 


E,(T) < Cl'(+o|| + (4Gc + C/-(r))||u-|b^(r) +4CcT). 


BP, 

P + 


Therefore, assuming that (|35l) is fulfilled and taking smaller c if needed, we get 


.n+l I 


Ep(T) < ■A+"'llRp(r) + C'^dkoll . g_i + 4CcT) 


BP, 

P,1 
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and thus, if 

(37) h^EAT) < 2C^'(lko|| . + 4 cCT) 

then also satisfies (f37|) . 

To complete the proof, we still have to justify that (j35h is fulfilled at rank n + 1. From 
the definition of || • \\ep{t)^ embedding and ([371) (at rank n + 1), we know that there exists 
some constant C" so that 


U^+\T) <C"(\\uo\\ d i+cT). 

K.i 

Hence there exists a constant c' > 0 such that if T and uq satisfy 
(38) ||mo|| d_i+cr<c' 

then both (1351) and ([37)1 are fulfilled at rank n + 1. 

If lliioll d_i > c' then we split vP into v^l+u^ with u^{t) Denoting ul '=^ e^'^uo 

and observing that uj = SnUi, we discover that 

U^{T) < \\Vul\\ d + ||VfI”|| d < C\\Vul\\ d + IIVS^II d . 

L}p{BlA 

The term with ul goes to 0 for T tending to 0 with a speed of convergence that may be 
described according to (fT^ . To handle the second term, we observe that satisfies 

= -S” • Vu" - ul • VS” - ul ■ Vul - V(G(a”)) - /(a”)^u”. 

Because S”'''^(0) = 0, combining (l30[) . product laws in Besov spaces and ((371) . we get 


(39) \\u^^^\\ep(t)<c( r \\u^ ._J|Vn”|| d dt + [^ 

\Jo Bf, Bl, Jo 


\Ul\\ d \\u 


dt 




rf_i T ^ ll^^ll ——1 T T||n^ 

Arguing by interpolation yields for any /? > 0, 




f 

Jo 


\ul\\ . d Ike'll ,ddt< /3||ui;,|| 


BP, 

pA 


BP, 

pA 


(s;,i ) L}piBP^ ) 


+Cfi-^\\UL\\ 


Br^Bpi ) B'A’i.Bpi ) 


T ^ p,l 


Besides, as||nL|| d_, ^CIIuoH d_,, Inequality (l3^ implies that 


~n+l 


\ep(t)<c{{U^{T) + 


d_l + /3 
RP 
°p,i 


I + ||®”|| -d ) 11 'W '"’11 (T) 


+ ||'Wo|| . d_i ||iil| 


BP, 

pA 


■ ^+1. 


+ (T + ||ul| 


• v + 1 


)l|a^ 


BliBP, ) ) L^iBP^ 


Choosing /3 = l/(4C'||uo|| .d_j), remembering (l37[) and ([3^ (taking c smaller if needed). 
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we conclude that there exists C" so that 


|~n+l II ^ ^ ||~n II 

\U \\Ep(T) < 2 11 “ llEpiT) 


+ C'c. 


Hence ||n"’||£;p(T) < ‘IC’c implies that ^ 2C""c, and thus ([37|) is fulfilled at 

rank n + 1 if c has been chosen small enough. 


Finally, let us notice that there exists some T > 0 so that 

(d l\.. ■ C 




“1“ ll'^oll . d_i 


P,1 


Hence (fT^ ensures that (jiO]) is satisfied for this choice of T. 


Step 3: Convergence in the case 1 < p < d and d > 3. Let da"- — aP and 5uP 


u 


n+i _ couple {daP^^, 5uP^^) satisfies 


(41) 


dtda3^+^ + • Vda^+^ = dF^, 

dtda'^+^ - AdF^^^ = X) ■=! 


with (5Ff • Va’"+\ <^ 2 ^ divu’"+\ (JFg” -(1 + a’")div5u’", 

dG2 • V(5u", dG2 -du^ ■ Vu”. 

Owing to the first equation, one can perform estimates for {daP,duP) only in a space with 
one less derivative, namely in C([0,T];i?p^ ) x Fp(T) with 

FpiT) C([0, T]; h|”') n L\0, T; h|i). 

Now, using the same type of computations as in Step 3, we get for all t G [0,T], 

||&”+Hi)ll 4-1 < ll&"+'(0)|| d_,+C ['\\W+^\\ d (ll&’^ll + dAdr 


p,i 


BP, 

pA 


BP, 

p,i 


BP, 

p,i 


BP, 

p,i 


+ [\l + C\\a^ d +C||a'^+i|| 4)||(it”|| d dr. 
Jo BF BF BF 


P,1 P,1 P,1 

Using the bounds provided by the previous step, we thus get, taking c smaller if needed. 


(42) ||&“+1|| <||4."+i(0)||,4_.+i||4."|l .4-. +2||&”IIf,«). 

As in the previous step, bounding in Fp{T) follows from (pO]) and product laws. 

However, as less regularity is available, one has to make the stronger assumption 

(43) d > 3 and 1 < p < d. 

Taking c smaller if needed, we eventually get thanks to (l34)l . (j35]) and (l37)l 

WSa-^^hpiT) < C||<5u"+'(0)|| + i(||&"|| + ||<5 u"||^,(t)). 

Bh 8 L^iBl, ) 
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Combining with (I42p yields 

ll&"+'ll +4 ||S<"+'IIf,(t) <C(||&”+‘(0)|| +4||S,"+1(0)|| 




^PT 






n I 


Ol o-riU O.’l' _ ,,0 


d_i +4||<^''||^^(7.)). 


Summing up over n E N, we conclude that (a” — a^)nGN and {u^ — u^)neN converge in 
. 4-1 

C{[0,T]-, Bp ^ ) and in Fp{T), respectively. 

Step 4- Checking that the limit is a solution and upgrading regularity. From Step 3, we 
know that there exists a and u so that 

. ^-1 

a"'— oo —>■ a — oo in L°°(0,T;i?p^ ) and u” — uq —)• u — uq in Fp{T). 

The bounds of Step 2 combined with Banach-Alaoglu theorem imply that in addition 


^ a 




in L°°{0,T-, Bp i) weak * and ^ u in L°°(0,r;i?p^ ) weak *. 
Routine verifications thus allow to pass to the limit in (|32p . 

d 2 

The previous step tells us that is in L^{0,T] B^ ^ ). To upgrade the regularity 
exponent by 1, let us write that for all J E N: 




u\\lp dt < 


2 ^ 1 - 1 )] 




j V u, IILP 

T 


dt 


+24 j:/ 


2^'^p"^^||AjV\ - AjV^u^Lp dt. 


The first term is bounded by the r.h.s. of p37p while, by virtue of Step 3, the second one 
tends to 0 when n goes 0. Hence, letting J tend to -|-oo ensures that ||V^u| 


is finite. 


Next, as (a, u) satisfies 


^H^pA ) 

~ . 4 

Theorem 12.21 and Inequality (1301) imply that a E L’^{Bpi) 


■ 


and that u E ^ ), which, combined with the fact that a — oq E C([0,T]; Bp^ ) and 

u — uo E C([0, T]; ) implies that a E C([0,T]; and u E C([0,T]; ). 

Step 5: Uniqueness. Consider two solutions (a^,u^) and of (IMp with the above 

regularity. The difference {6a, du) {a^ — a^,u^ — u^) satisfies 

dtda + ■ V6a = 

dtdu - A6u = Yl\=i 


(44) 


with 6 Fi = —Su-Va^, 6 F 2 = —5idivu^, dFs = —(1-|-a^)div&t, 

6 Gi (/(a^) - I{a'^))Au‘^, 6 G 2 -I{a^) A 6 u, 6 G 3 V(G(ai) - G{a‘^)), 

dGi -u^ • V 6 u, 6 G 5 -(5u • Vu^. 

Mimicking the computations of Step 3, it is easy to see that if (l43P is fulfilled then { 6 a, 6 u) = 
0inC([0,r];H|“')xFp(r). 
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It turns out that the limit case d = 2 or p = d tractable even though 
(45) G i and Au^ G ^d,i} implies SaAu^ G only. 

Now, applying Theorem 12.21 and product laws (see Proposition 12.3p gives 


+ (1 + IIo IIl?? 




Regarding du, owing to (I45p . one has to apply Remark 12.21 rather than Theorem l2.ll which 
enables us to control the following quantity: 


ll^ll 


def 




= sup 2-^ IIA 




which is slightly weaker than ||(5u||^i 


Inserting the following logarithmic interpolation inequality (see my- 


(46) ||<iu||^^(^i^^) < ||&x||j^(^i^^)log(^e + 


11*^11 rl 


'^tKoo) 

in the estimate for 6 a and using Osgood lemma (see e.g. [T], Chap. 3), we end up with 

,exp(-/o Odr) 


ll&ll 




+IIM 




m +||(it(o)||^-i 


where a is in T^(0, T) and depends only on the high norms of the two solutions. This yields 
uniqueness on [0,T]. 


3.2. The Lagrangian approach. We now propose another proof of the local well-posedness 
of (|26h , which will provide us with the statement of Theorem 13.11 in its full generality. It is 
based on the Lagrangian formulation of the system under consideration. 

To make it more precise, we need to introduce more notation. First, we agree that for a 
function F : ^ x R”* then divT : R'^ —)• R”^ is dehned by 

d 

{divFy ^ diFij for 1 < j < m, 

i=l 

and that for A = (^ij)i<i,j<d and B = {Bij)i<ij<d two d x d matrices, 

A : B = Tr AB = AijBji. 

ij 

The notation adj (^) designates the adjugate matrix that is the transposed cofactor matrix 
of A. Of course if A is invertible then we have adj (A) = (det A) A~^. Finally, given some 
matrix A, we define the “twisted” deformation tensor and divergence operator (acting on 
vector fields z) by the formulae 

Da{z) - {^Dz ■ A+ '^A - Vz) and div^ z'^='^A-. Vz = Dz : A. 

We recall the following classical result (see the proof in e.g. m)- 
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Lemma 3.1. Let K be a sealar function overW^ and H, a veetor-field. Let X be a 
diffeomorphism such that J ='^ dei^DyX) > 0. Then we have 

(47) VJL = J-U\Yy{v.A^{DyX)K), 

(48) = J-^ divy (adj {DyX)H). 

Let X be the flow associated to the vector-field u, that is the solution to 

(49) 

def 

Let Q{t,y) = g{t, X{t,y)) and u{t,y) = u{t, X{t,y)). Formally, we see from the chain rule 
and Lemma [3lT] above that (^?, u) satisfies (jH) if and only if {g,u) fulfills 

dt{Jg) = 0 

— div ^adj (L>X) (2^(^)L)yi(M)-|-A(^) divA uid -|-P(^)Id)^ =0 


X{t,y)=y+[ u{T,X{T,y))dT. 
Jo 


(50) 


def 


def 


with J det DX, A (DyX) ^ and 

(51) X { t , y ) = y+f u { T , y ) dT . 

Jo 

The first equation means that g = and the velocity equation thus recasts in: 

Lgoiu) = ^>0 Miv(/i(tt,tt) -h h{u,u) -I- h{u,u) + h{u)) 

with 

(52) Lyg (u) dtu - £>0 Miv {2g,{go)D{u) + A(£»o)divu Id) 
and 

Li{v,w) (adj (DXy) - ld){fi{JT^go){DwA, + ^A,Xw) + A(JT^^?o)(^A : Vu;)Id) 

h{v,w) go) - yi{gQ)){Dw A^ + '^A^Vw) + go) - \{gQ))(^A^ : Vr(;)Id 

h{v,w) /i(£io)(Zltc(A - Id) + ^(A - Id)Vt(;)-b A(^>o)(’^(^t; - Id) : Vr(;)Id 
Li{v) -adj {DX,)P{goJT^), 

where Xy is given by ([4^ with v instead of u, Ay {DXy)~^ and Jy dei DXy. 

So finally, in order to solve (f50]l locally, it suffices to show that the map 

(53) <h : V I—^ u 
with u the solution to 

7^eo(^) = £’o^div(/i(u,u) + L2{v,v) + L3{v,v) + L^iv)), 

n|t=o = no 

has a fixed point in Ep{T) for small enough T. 

In order to treat the case where g is just bounded away from zero, we need to generalize 
(j30p to the following Lame system with nonconstant coefficients: 

(54) dfU — 2adiv {fj,D{u)) — 6V(Adivn) = /, 
where a, 6, A and y satisfy the following uniform ellipticity condition: 

=^min(" inf (ou)(t,x), inf (2aa + bX)(t,x)) >0. 
V(La:)G[0,r]xB‘^ xR'^ / 


(55) 


a 


(t,x)g[0,T] xl 
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In [16], the following statement has been proved. 

Proposition 3.1. Let a, b, A and fa he bounded functions satisfying (|55p . Assume that 


■ 


aV/i, 6VA, gS/a and Wh are in L°°{0,T] ) for some 1 < p < oo. There exist r] > 0 

and a > 0 such that if for some m £ Z we have 


(56) 


mm 


( inf 5m(2a/i + 5A)(t,x), inf Sjn{ap){t,x)j > — 


V(t,x)g[0,T]xK‘^ (La:)G[0,r]x 

(57) ||(Id -5m)(/iVo,oV^,AV6,6VA)|| d_,<ria 

) 

then the solutions to (l5^ satisfy for all t G [0,T], 

+ Oi\\u\\Ll(B;ff) 

rt 


\u\ 


LriBti) 


<C{\\uo\ 


B: 


pA 


ri(^s -,)expr— [ ||5m(/xVa, aV;U, AV6, 6VA)||^ d dr 

V«io 

P,1 


whenever — min(d/p, d/p') < s < d/p — 1. 

In order to show that $ in (|53l) admits a fixed point in Ep(T), we introduce, as in the 
previous subsection, the solution ul in Ep{T) to 

Liul = 0, u\t=o = Uq. 

We want to apply Banach fixed point theorem to in some suitable closed ball Bep(t){ul, R). 

def ~ def 

Let V be in Bep(t){'^LiR) and u = <h(i;). Denoting u = u — ul, we see that 


(58) 


LqoU = Qq Miv(/i(i;,i;) + hiv^v) + h{v,v) + h{v)) + (Li - Lg^)uL, 

u\t=Q = 0. 


The existence of some m G Z so that 
min I inf Sm 


n f^(go) ^(go) 


go 


go 


, inf 


a 

>2 


and 


(Id - Sm) 


(P{qo)^^ /^'(go)v7„ '^(go)v7„ 

^ 9 ^^0? ^^0? 9 ^1^0 


go 


go 


go 


7 /^(go) \ 
^ go ^ 


go 


d_, < rja 

ttp 

^p.i 


is ensured by the fact that all the coefficients (minus some constant) belong to the space 

. i 

Bp I which is defined in terms of a convergent series and embeds continuously in the set of 
bounded continuous functions. Hence, if one can show that the right-hand side of (1581) is 
. --1 

in L^(0,T;Hp^ ) (which will be carried out in the next step) then we will be allowed to 
apply Proposition 13.11 to bound u in Ep{T). 

First step: Stability of Bep{t){ul, R)- Let v G Bep(t){ul, R) and u be given by (j58ll . Let 
«o = go ~ 1- Proposition 13.11 product laws in Besov spaces and Proposition 12.41 imply that 


(59) ||tt||£;p(T) < Ce^‘^o'^{ ||(Li - Lp^)ul\\ 


LfiBE ) 


+ (l + Ikoll . d ) (||.f^4(?^)|| .d ^ ||/j(u, u)|| d )') 


for some constant depending only on go- 
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In what follows, we assume that T and R have been chosen so that, for a small enough 
positive constant c, 

rT 


(60) 


/ llVHI.i 
Jo 


dt < c. 


Now, using the decomposition 

(Li - Leo)uL = - l)div {2 ^{qq)D{ul) + A(£»o)divttL Id) 

+div(2(/i(£io) - ^^{l))D{u) + (A(£>o) - A(l))divitld), 


■ 


and Proposition 12.41 we see that {Li — Lg^)uL G L^(0,T;i?p^ ) and 


(61) 


||(Li-L,Jui||. 


. -2^-1 ^ 


IIroII d (1 + Ikoll d )II-C>'*^lII d 

1 (nP 


Likewise, flow and composition estimates (see the appendix) ensure that 

(62) ||/i(i;,'u;)|| ^(l + ll«o|| d )\\Dv\\ a \\Dw\\ d fori = 1,2,3, 

(63) ||/4(i;)|| , .. <T(l + ||ao||., )(l + ||Ilii|| , .. ). 

So plugging the above inequalities in ([59]) and keeping in mind that v satisfies ([60|1 . we get 
after decomposing v into v + ul- 


u\\ep(t) ^ Ce^‘^0’^'-^ (1 + ||ao|| , d 1{T + 


BP, 

Pil 




) 




+ \\DuLf d +{\\Dul\\ d +\\Dv\\ d )\\Dm d 

LUBl,) L^B^,) L},{Bl^) L},{Bl^) 


Now, because v € BEp{T)id, R), 

\\u 


E (T)<Ce^Po'^{l + \\aQ\\ d )^((r + ||oo|| d \\Dul\\ d ) 

BP- V RG L\.{Bl^) 


P,1 


+ (P + IIL^HlII .d )||Z)tti|| d -\-R^ 


LUKi) 

Therefore, if we first choose R so that for a small enough constant r] 

(64) (1 + llaoll dfR<r] 

and then take T so that 

(65) Cg^T <\og2, T<R^, ||ao|| d ||L>?xl|| d <R^, \\Dul\\ 


LHBP,,) 


Lh(Bli) 


d < P, 


then we may conclude that 4> maps Bep(t){'^l, R) into itself. 


Second step: contraction estimates. Let us now establish that, under Condition (1651) . the 
map d* is contractive. We consider two vector-fields and in Bep{t){'^l-,R), and set 

y} $(u^) and v? <I>(u^). Let 5u‘^= v? — v} and 5v‘^= — v^. We have 

Lg^Sa = ^>0 Miv(^(/i(i;^,i;^) - /i(i;\u^)) 

+{h{v\v‘^) - h{v\v^)) + {h{v\v^) - h{v\v^)) + {h{v^) - h{v^)))- 
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So applying Proposition 13.11 frecall that < log 2), we get 


3 

(66) \\Su\\ep(t) < (^(1 + ||ao|| . | ) 

i=i 




+ WUiv^) - Uiv^] 




In order to deal with the first term of the right-hand side, we use the decomposition 

Ii{v‘^,v‘^) - = X{J~2^go)(^A^2 :Vv‘^){adj{DX^2) - adi{DX^i)) 

+ (adj {DX,i) - Id) (A(J;,ipo) - A(J-i^po)) i^A ,2 : Vv^) 

+ (adj {DX,i ) - Id) A( go) {CA ,2 - ^A,i ) : + ^^,2 : X6v) 

+ terms pertaining to fi. 


Taking advantage of product laws in Besov spaces, of Proposition 12.41 and of the flow 
estimates in the appendix, we deduce that for some constant Cg^ depending only on go: 

,d <Cg,\\iDv\Dv^)\\ ,d IID&II . 

Similar estimates may be proved for the next two terms of the right-hand side of (I66h . 
Concerning the last one, we use the decomposition 

h{v^) - h{v^) = (adj {DX,i)- adj (I)X, 2 ))P(J-/po) - adj {DXpi){P{J; 2 ^go) - PiJ^i^go))- 
Hence 

Whiv^) - hiv^)\\ d <C(l + ||ao|| . )r||I)<fc|| .. 

BP, 

We end up with 

\\du\\EpiT)<C{l + \\ao\\,d)^{T+\\iDv\Dv^)\\ )\\D6v\\ . 

B^i LUB^i) L}p(Bl,) 

Given that and are in Bep{t){ul, R), our hypotheses over T and R (with smaller rj 
in (IMl) if need be) thus ensure that, 

\\Su,\\ep(T) < -^W^WEpiT}- 

Hence admits a unique fixed point in Bep{t){ul-,R)- 


Third step: Regularity of the density. Set g By construction {g^u) satisfies (IHH]) 

and a '=^ p — 1 is given by 

a = {Ju ^ ~ l)®o + ao- 

_ d _ d 

From the appendix, as Du G L^(0,T;Hp^), we have — 1 belongs to C([0, T]; ^). 

d d 

Hence a is in C([0,T]; Bpi), too. Because B^^ is continuously embedded in L°°, the density 
remains bounded away from 0 on [0,T] (taking T smaller if needed). 
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Last step: Uniqueness and continuity of the flow map. Let the data {qq,Uq) and 
fulfill the assumptions of Theorem EH and let and be the corresponding 

solutions. Setting 5a^ , we see that 


LsliSu) = (L^i -L^ 2 )(n 2 ) + (qI) Miv n^) -/|(n\ + (/|(u2) -/|(u^))) 

i=i 

+{Ql)~^diY - lj){u^,u^) + (/| - /l)(n^)), 


i=i 


where I\, /| and II correspond to the quantities that have been defined just above 

(El, with density Qq. Compared to the second step, the only definitely new terms are 
{L^i — L^ 2 ){u^) and the last line. As regards (L^i — L^ 2 ){u‘^), we have for t <T, 



The other new terms satisfy analogous estimates. Hence, applying Proposition 13.11 yields 
if 6po is small enough: 

\\du\\E^(t) <C^i{{t+\\Du^\\ d +\\HEp(t))\M\Ej,(t) 

+ ||&io|| d + ll'^Poll d (^ + ll-^^^ll d )). 

An obvious bootstrap argument thus shows that if t, Suq and 5po are small enough then 

||-5u||Ep(t) < 2C'go(||(5uo|| + ll^/Ooll )■ 

1 1 

As regards the density, we have da = + {J~ 2 ^ — J~i^)aQ. Hence for all t G [0,T], 

||<fa(t)|| d <C(1 + ||W|| d )\\dao\\ d\\Ddu\\ d . 

So we get uniqueness and continuity of the flow map on a small time interval. Then iterating 
the proof yields uniqueness on the initial time interval [0, T], as well as Lipschitz continuity 
of the flow map. 

It is now easy to conclude to Theorem 13.11 in its full generality, as a mere corollary 
of the following proposition which states the equivalence of Systems (f26]l and ([50l) in our 
functional setting (see the proof in US]). 

. ^ 

Proposition 3.2. Let 1 < p < 2d. Assume that the couple {g, u) with (^>— 1) G C([0, T]; 
and u G Ep{T) is a solution to (IM]l such that 

(67) [ ||Vrt|| d dt < c. 

Jo Bfi 

j— ]— i 

Let X be the flow of u defined in dH. Then the couple {g, u) = {g o X,uo X) belongs to 
the same functional space as {g,u), and satisfies (l50l) . 
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Conversely, if {q — l,u) belongs to C{[0,T]] B^-^) x Ep(T) and {q,u) satisfies ([50l) and, 
for a small enough constant c, 

(68) [ llVtill d dt < c 

Jo Bf, 

then the map Xt ='^ X{t, •) defined in d^T]) is a diffeomorphism on and the couple 
{g,u){t) =‘^ (git) o Xf^ ,u{t) o Xfi^) satisfies (I^SD and has the same regularity as {g,u). 


4. The global existence issue 

This section is devoted to the proof of global existence of strong solutions for small 
perturbations of the constant state (g, u) = (1,0), under the stability assumption P'{1) > 0. 
For simplicity, we assume that the viscosity functions A and p, are constant. 

Let us emphasize that the approach we used so far to solve (f2U|l cannot provide us with 
global-in-time estimates (even if both oq and uq are small) because we completely ignored 
the coupling between the mass and momentum equation through the pressure term and 
looked at it as a low order source term, just writing 

dtu — Au = —u ■ Xu — I{a)Au — V(G'(o)). 

Then, applying Inequality (f30]l and product laws in Besov spaces led to 

(69) < C'f Ikoll . d_i + ||u||e (t)(||a|| d + \\u\\E„{t)) + \\a\\,ddT 

V LrlBl,) Jo Bf, 

At the same time, as a is a solution to a transport equation, we can only get bounds on 
||a|| d and the last term of (l69l) is thus out of control for t —>■ -|-oo. 


4.1. The linearized compressible Navier-Stokes system, and main result. The key 

to proving global results is a refined analysis of the linearized system (IM]l about (a, u) = 
(0,0) taking the coupling between the mass and momentum equation through the pressure 
term into account. The system in question reads: 


(70) 


f dta + divM = /, 

\ dtu — pAu — (A -|- p)Xdivu + P'{l)Xa = g. 


Applying the orthogonal projectors V and Q over divergence-free and potential vector- 

fields, respectively, to the second equation, and setting a P'{^) and u A -|- 2p, System 
(EQD translates into 


I dta + dwQu = f, 

(71) < dtQu - uAQu + aXa = Qg, 

I dtVu — pAVu = Vg. 


We see that Vu satisfies an ordinary heat equation, which is uncoupled from a and Qu. 
For studying the coupling between a and Qu, it is convenient to set v |iA|“^divtt (with 

T'(|Z1|^m)(^) |?|'^S(0)) keeping in mind that, according to (fTTIl . bounding v or Qu is 

equivalent, as one can go from v to Qu or from Qu to v by means of a 0 order homogeneous 
Fourier multiplier. 
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For Rotational simplicity, we assume from now on thalQ a = V = 1. Hence (a, u) satisfies 
the following 2x2 system: 


(73) 


dta + \D\v = /, 

dtv — Av — \D\a = /i |i7|“^divg. 


Taking the Fourier transform with respect to x, and denoting p |^| with ^ G the 
Fourier variable, System (f73ll translates into 


(74) 



with Mp 


0 -p 
p 


• In the low frequency regime p < 2, Mp has two complex conjugated eigenvalues: 

A±(p) - y(1 ± iS{p)) with S{p) - 1 

which have real part —p^/2, exactly as for the heat equation with diffusion 1/2. 

• In the high frequency regime p > 2, there are two distinct real eigenvalues: 

A±(p) =^-^(1 ±72(p)) with R{p) - A . 

As 1 — ^ ^jfP' for p — y -|-cx), wg hsiVG A_|_(/^) ~ —p^ and X-{p) — 1- In other 

words, a parabolic and a damped mode coexist. 

Optimal a priori estimates may be easily derived by computing the explicit solution of 
(j73l) explicitly in the Fourier space. Below, we present an alternative method which is 
generalizable to much more complicated systems where explicit computations are no longer 
possible (see e.g. m)- 

Fix some p > 0 and consider the corresponding solution {A,V) of (1741) in the case 
f = h = 0. We easily get the following three identities: 

(75) l^|yi|2 + pRe(HF) = 0, 

(76) i||y|2 + Y|Fp-pRe(HF) = 0, 

(77) ^ Re (HR) + p\V\^ - p\A\^ + p^ Re (HR) = 0, 
from which we deduce 

(78) \j^Cl + p\A,V)\^ = 0 with £2 drf2|(H,R)|2 + \pA\^ - 2pRe(HR). 

Using Young inequality, we discover that there exists some constant Cq > 0 independent 
of p so that 

(79) C^^Cl<\{A,pA,V)\^ KCoCl 

Combining with (|78p . we conclude that there exists a universal constant cq > 0 so that 

(80) < g-comin(i,p 2 )t^ 2 (Q^ forall t>0. 


'^Which is not restrictive as the rescaling 
(72) a{t,x) = t, ^ x) and u{t,x) = t, ^ x) 

ensures that (a, u) satisfies dZH) with a = V = 1. 
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In the case of general source terms / and h in (j74j) . using Duhamel formula and applying 
the above inequality to A{t) = a{t,p) and V{t) = v{t,p) thus leads to 

(81) \{a,pa,v){t)\p‘^)J \{a, pa,v)\dr < C (^{ao, pao,vo)\ + \{f,pf,h)\dT 

Note that as 


dtv + p‘^v = h + pd, 


we also have 


p^ f |u(r)| dr < |uo| + [ |/i(r)|(ir+ f |/9a(r)|dr. 
Jo Jo Jo 


And thus, bounding the last term according to (1811) . we get the following inequality which 
provides the full parabolic smoothing for v: 


(82) |(d, pd, u)(f)| + min(p, p^) f |d| dr + f |d| dr 

Jo Jo 

< C'^|(do,pdo,do)| + ^ |(/,p/,/i)| dr 

From Inequality (|80p and Fourier-Plancherel theorem, it is easy to obtain estimates of 
type for the solutions to (j70p . Optimal informations will be obtained if splitting the 
unknowns into frequency packets of comparable sizes. To this end, one may apply to 
(f70|) and get 

( cljAfca + divAfcQu = Afc/, 

(83) < dtAkQu-AAkQu + VAka = AkQg, 

[ dtAkVu - pAAkPu = AkPg. 

In the case with no source term then using (|80l) combined with Fourier-Plancherel theorem 
readily yields for some universal constant Cq, and cq depending only on p, 

||(Afca,AfcVa, Afen)(t)|A2 < AfcVa,Afcn)(0)|A2. 

Then, for general source terms, using Duhamel’s formula and repeating the computations 
leading to (l8^ . we end up with 


dr 


( 84 ) ||(Afca, AfeVa, Afcu)(t)||£,2 -Fmin(l,2^) /" || AfeVo||£,2 dr2^^ / ||Afctt|| 

Jo Jo 

< C^IKAfeOo, AfcVoo, AfcUo)||L2 -h j ||(Afc/, AfcV/, Afcp )||£,2 dr 


Multiplying both sides by 2^^, taking the supremum on [0,t] then summing up on A; > ko 
or k < ko, we conclude to the following: 


Proposition 4.1. Let s G M and {a,u) satisfy (I70p with P'{1) = v = 1. Let ko G Z. Then 
we have for some constant C depending only on ko and p, and all t >0, 

\\i 


[a,u 


Ilf 




+ 




2,1 


<C(||(a„.a„)||^ +||(/,9)||'^,„ ), 
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where we used the notation 


(85) 



2 ^-\\Kkz\\Lr and 2'^"||Afcz||LP. 

k<ko k>ko 


The high frequencies inequality means that in order to get optimal estimates, it is suitable 
to work with the same regularity for Va and u. In contrast, for low frequencies, one has to 
work in the same space for a and u, a fact which does not follow from our rough scaling 
considerations ([5]) but is fundamental to keep the pressure term under control in (I26p . 

Granted with the above proposition, it is now natural to look at (f26ll as System (TTn]) 
with right-hand side 

/ = —div(au) and g =—u-Vu — I{a)Au — k{a)Va where k{a) G'{a) — G'{0). 


The problem is that / will cause a loss of one derivative as there is no smoothing effect for 
a in high frequency. A second limitation of Proposition 14. II is that it concerns Besov spaces 
related to whereas we know the system to be locally well-posed in more general Besov 
spaces (see Theorem ET]). To overcome the hrst problem, let us include the convection 
terms in our linear analysis, thus considering: 


( 86 ) 


J dta + V ■ Va + divu = f, 

\ dtu + V ■ Vu — gAu — (A -|- g,)Vdivu + aVa = g, 


where v stands for a given time-dependent vector field. 


Proposition 4.2. Let —d/2 < s < d/2 and {a,u) satisfy ((86]) with a = u = 1. Let k^ E 
Then we have for some constant G depending only on /cq and g, and all t > 0, 

I k 


(“, + 11 11G (RIP + II“IIg(R|+P 


<C( \\iao,Vao,uo)\\jSs+\\if,Vf, 


2,1 


IG(B|,i) 


+ / l|Vu|| d ||(a, Va,u)||o. dr 


^ 2,1 


2,1 


Proof. Applying to (f86]l yields: 

f dtOkP Ak{v-Va) + divUk = fk, 

\ dtUk + Ak{v ■ Vu) - gAuk - (A -F g)Vdivuk + Vok = gk, 


with Ok = Afco, Uk = AkU, fk = Akf and gk = Akg. 

Keeping in mind the proof of Proposition 14.11 we introduce 

— 2|| (Ofc, I^A:)|li2 + ||Vafc||^2 + {uk I Vak)L2. 

Now, remembering that A -|- 2|U = 1, we get 

(87) T\\'^'Puk\\‘^j^2 + ||(V Qufc, Vafc)||^2 = {gk \ {2uk+Vak)) j^2 + 2(/fc | ak)L2 

+ (V/fc I Vafc)^2 - 2(Afc(u • Va) \ Ok)^2 - ‘2^{^k{v ■ Vn)|nfc)^2 

{AkV{v ■ Va) I Vak)j^2 - (Afc(u • Vu) \ Vak)^2 - (AfcV(u • Vo) | Uk)^2- 

Let us explain how to bound the convection terms. To handle the second and third terms 
of the second line, we proceed as explained below, taking b E {a, u^, - ■ ■ , 
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Integrating by parts and setting 


def 


Rk{v,b) = Ak{v-Vb)-vVAkb, 

we discover that 

{Ak{v-Vb) \bk)j^2 = j{v-Vbk) ■ bkdx + J Rk{v,b)bkdx 

- j l^kl'^divvdx + \\Rk{v,b)\\i2\\bk\\L^. 

Bounding the last term according to (I24h . we thus get 


|(Afc(u-V6) I bk)^2\<Cck2-'^^\\Vv\\ , I 

^21 




Whh^ 


def 


with {ck)k& in the unit sphere of £^(Z). 

Next, we use the fact that for i G {1, • • • , d}, 

diAk{v ■ Va) = V ■ Vdiak + Rk{v,a) with Rl.{v,a) [diAk,v] ■ Va. 

By adapting the proof of (l2^ . it is easy to prove that 

\\Rk{v,a)\y<Cck2-’^^\\Vv\\ , ||Va||g 

H ^ 2,1 

^2,1 

Then using an integration by parts, exactly as above, we conclude that 

{AkV{vVa)\Vak )^2 < Cck2-^^\\Vv\\ . , ||Va||^. J|Vafc||i 2 . 

^2,1 

Finally, to handle the last two convection terms, we use the fact that 

{Ak{v ■ Vu) I Vak)^2 + (Aa:V(u • Va) | Uk) ^2 

= {v ■ Vufc|Vafe )^2 + {{v ■ V)Vafc)|ufc )^2 + {Rkiv,u) \ Vafc )^2 + {Rk{v,a) \ uu)^ 2 - 

Integrating by parts in the first two terms of the second line and using again (|24p to bound 
the last two terms eventually leads to 

[Ak{v • Vu) I Vafc )^2 + (AfcV(u • Va) | Uk )^2 

< Cck2-^^\\Vv\\ . . (IIVail^ ||ufc||^2 + ||u||^. ^ ||Vafc||i2). 
^ 2,1 

Because Ck ~ IK^fc) Vafc,Ufc)||^ 2 , we thus conclude that 

(88) + ^l\\VVuk\\l2 + IKVQufc, Vafc)||i. 

< (ll(/fc,V/fc,5fc)||L2+Ccfc2-"^||Vu|| . ||(a,Va,u)||5. )£fc, 

^ r>2 ^2,1^ 

^2,1 

which after time integration and multiplication by 2^^ yields 

2'^^£fe(t) +co2'^^min(l,22^) T ||(a^, Va^, Ufe)||i 2 dr < 2^^Lk{£)) + f 2 ’^^\\gk\\L 2 dr 

do do 
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Taking the supremum on [0, t] then summing up over k, we thus get 


(89) + \\{a,u)\f^s +2 dr + ||(Va, u)||^, ^ dr 

<||(a,Va,u)(0)||^. + f IK/, V/, 9)11^. dr + [ ||Vu|| d ||(a, Va,u)||g. dr. 

Jo Jo ^ 2.1 ’ 

Finally, using the fact that 

dtu + V ■ Vu — fiAu — (A + ^)Vdivu = g — Va, 
localizing according to A^, and arguing as above, we find out that 




f 




|u(0)|| 


BS 


+ [ lid-Va||^.^dr+ f ||Vu|| d ||u||^s^dr. 
Jo ' Jo ^ 2,1 ’ 


Then bounding Va according to (j89l) completes the proof of the proposition. □ □ 


It turns out to be possible to extend the above proposition to more general Besov spaces 
related to the spaces with p ^ 2. The proof relies on a paralinearized version of System 
(I86p combined with a Lagrangian change of variables (see HE!). Here, in order to solve 
(j26p globally, we shall follow a more elementary approach based on the paper by B. Haspot 
m : we use Proposition 14.11 only for bounding low frequencies, and perform a suitable 
quasi-diagonalization of the system to handle high frequencies. This eventually leads to 
the following statemenl@ that will be proved in the rest of this section: 


Theorem 4.1. Let d > 2. Let p G [2, min(4, 2d/(d — 2)] with, additionally, p 7 ^ 4 if d = 2. 
Assume with no loss of generality that P'{1) = 1 and n = 1. There exists a universal integer 

— . — — 1 

feo £ N o-nd a small constant c = c{p,d, p,G) such that if oq G Bpi and uq G with 

■ --1 

besides (qqjUq) in Bf^ (with the notation = Sk(,+iz and = z — z^) satisfy 


(90) 


=^||(ao,^o)||^ , + l|aofd +||nofd ,<c 
^21 


then ([ 26 ]) has a unique global-in-time solution (a, u) in the space Xp defined by 


1 


{a,uy G a(M+;B|, ) n T'(M+;B|K ), G CyR+-,By,) n L^(M+;B;^i), 




eCb{R+;B;, )nL^(M+;B;, ) 


where we agree that Cb{R^; B^ y B^ y n L°°{M.j_; B(i y, s G M, 1 < g < 00. 
Furthermore, we have for some constant C = C{p, d, p, G), 

( 91 ) ||(o,a)||x, <C'Xp,o. 


Remark 4.1. Condition (j90p is satisfied for small oq and large highly oscillating velocities: 
take Uq : X i-A (f>{x) sin(e“^x • oj) n with u and n in and (j) G 5(M'^). Then 

\\uo\\ <Ge^~p ifp>d, 


’The reader may refer to [18] for a slightly more general result. 
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and llugll^ c 

b:^ 


has fast decay with respect to e. Hence such data with small enough e generate 


global unique solutions in dimension d = 2,3. 

Remark 4.2. One may extend the above global result to 2d/{d + 2) <p <2 provided the 
following smallness condition is fulfilled: 






d + lkolL#_i < V- 


B.i 


'' 2,1 ' '-^ 2,1 ^ 2,1 

Indeed, Theorem \5.1\ vrovides a global small solution in X 2 . Therefore it is only a matter of 
checking that the constructed solution has additional regularity Xp. This may be aehieved 
by following Steps 3 and 4 of the proof below, knowing already that the solution is in X 2 . 
The condition that 2d/{d + 2) < p comes from the part ■ Va of the convection term in 

. d 

the mass equation, as is only in Bf ^), and the regularity to he transported is 

. i 

Bp^. Hence we need to have d/p < d/2 + 1 (see Theorem \2.Sfl . The same condition appears 
when handling k{a)Va. 

Remark 4.3. Using spaee C;)(]R+;-B| ;^) rather than just Cb(^+', B 2 i) is not essential in 
the proof of Theorem \4.1\ We chose to present that slightly more accurate result, as it will 
be needed when investigating time decay estimates, at the end of the survey. 

4.2. Global a priori estimates. Consider a smooth solution (a, u) to (|26p satisfying, say, 

(92) ll«llL°°(R+xRrf) < 1/2- 

We want to find conditions under which the following quantity: 


\\{a,u)\\i d + 
LTiBh ) 


[a,u) 


d-Ul 


+ ki 


LriBfi) 


+ iior +11^^11: 


. ^ — 1 

'(Ki ) 


+ rw 


+ 1 . 


LUBf, ) 


satisfies (I9ip for all t G 
Rewriting System 


as follows: 


f dta + divu = f—div{au), 

\ dtu — Au + Va = 5 '=^ —u ■ Vu — I{a)Au — k{a)Va, 

we shall take advantage of Proposition 14.11 with s' = d/2 — 1 to bound the low frequency 
part of {a,u). To handle high frequencies, following [22], we shall use the facts that, up to 
low order terms: 

• Vu satishes a heat equation (hence parabolic smoothing in any Besov space); 

• The effective velocity 


(93) w = V(—A) ^(o —divu) 

satisfies a heat equation; 

• The high frequencies of a have exponential decay. 


First step: Low frequencies. From Proposition 14.11 we readily infer that 

(94) ||(a,u)||l .d , + 




ia,u)\\\ <||(ao,uo)||%_, + ||(/, 
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Second step: high frequencies, the incompressible part of the velocity. To handle Vu, we 
just use the fact that 

dtPu — ji/Wu = Vg. 

Hence, according to Remark 12.21 (restricted to high frequencies) 

(95) \\Vn\t +(<l|Pu|l'' <C(||Pu„||'‘,_, + ||P 9 ||'‘ 

) '-■(B'l ) B'l 1,>(B;, ) 


Third step: high frequencies, the effective velocity and the density. On the one hand, the 
effective velocity w defined in ([93l) fulfills 

dtw — Aw = V(—A)“^(/ — div( 7 ) +w — (—A)“^Va. 

Therefore, Theorem O and the fact that V(—A) ^ is an homogeneous Fourier multiplier 
of degree —1 imply that 


II lift def II Ii/i 

(96) \\wr~ = \\w\\f 




) 

On the other hand, we have 

(97) SfO + div(ari) + a = —divrc. 

We claim that 


+ ||/-div5f +ll«^-(-A)-'Vaf 




(98) ||a||^ d + l|a||^ i <C'(||ao||^d + ||divr(;||'' d + / UVuU d ||a||_d dr 


BP, 

pA 


LUBfi) Jo 

Indeed, as in the proof of Theorem 12.21 let us apply A^ to (fWD . We get 
dtA^a + u ■ VAfcO + A^a = — Afc(adivtt) — A^divtc + i?fc, 
where, according to ()24ll . the remainder term Rk satisfies: 

\\Rk\\Lp < Cck2~^p\\'Vu\\ d ||a|| d with = 1 

rS P rS P / ^ 


BP, BP, 

p,l p,l 


BP, BP, 

P,1 Pil 


k£Z 


and 


|adivw|| d < Clldivull d ||a|| d . 

^ - op 

pa 


^pa 


BP, 

pa 


Therefore evaluating the norm of A^a seen as the solution to a transport equation, 

_ 

multiplying by 2 p and summing up over k > ko yields (|98l) . 

Next, let us observe that, owing to the high frequency cut-off, we have for some universal 
constant C, 


\w\ 


d+i and ||(-A)-iVaf d_^ <C2 -2^°||af d . 


BP, 

pa 


-+i 

BP, 

pa 


BP, 

pa 


BP, 

pa 


(99) 
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In consequence, combining (IMIl and (1981) . and choosing ko large enough yields 


( 100 ) ll'R’II S + ||o(^) 11^ d “1“ ll®ll^ d < 

Bil LUbI,) 


+ 11 /- 


^ I lu 11 ^ 

d_i + ||fl0|| d 

■ p J - 

V 


r> P R P 

1 1 


divfifll^ d_2 + [ ||Vn|| d 
LliB^a ) Jo Bl, 


all d dr 

J^P 

p,i 


Fourth step: end of the proof of the linear estimate. Putting Inequality (llOOl) together with 
(fMl) and (1951) and observing that 


I“IIe,(i) s IIP” 


I ^ 

I Ep (t) 


+ 11 11 ? / j -\ + ( 11 ^ 11 ^ d_ 


I ~ 

Ep{t) 


+ llttl 




)> 


we come to the conclusion (if ko has been taken large enough) that 

Xp{t) < Ck, (xp{0) + /'(IK/, g)\f, , ^ + ll/f + llfff d_, + llVull . d ||a|| . d ) 

V Jo Bl, Bl, BP, BpP, BpP, J 

Fifth step : nonlinear estimates. It is only a matter of proving that under hypothesis ()92p . 
we have 

(101) [\\\if,g)\f. d , + ll/f d_, + Ibf d_, + llVull . d ||a|| . d ) dr < CX^it). 

J^ Bl Bf, Bf, S/i BpPi 

As p > 2, it is clear that the last term of the r.h.s. of (I101|) is bounded by CXp{t). Next, 
arguing exactly as in the proof of the local existence, we easily get for 1 < p < 2d, 

ll/ll d_, <C'||a|| d ||m|| d , 

LHbE ) LUBf,) L?(Bpf) 

f II . d_, < Cdlull . d_, ||Vn|| d 

LliBli ) ) LliBf,) 

+ lfII — 11^ ^11 — — 1 + IfII — ll^^ll — — 1 )* 

LFiBfi) LliBf^ ) Lf(BP,) LUBf, ) 

Therefore, using the definition of Xp(t) and embedding (recall that p >2), we get 

ll(/,ff, <CX^p{t). 


LliBL ) 


So we are left with the proof of 
( 102 ) 


IK/,fir .d, <cx^{t). 

f (^2^ ) 

Let us admit the following two inequalities (the first one being proved in [18] and the second 
one being a particular case of Proposition 12.31 followed by suitable embedding, owing to 

1 <p/2 < 2) : 

(103) llTaftll .,_i+d_| < (Pllall |_,||6||^.^ if d>2 and ^ < p < min(4, , 

^2,1 Bpi P’ 

(104) ||i?(a,6)|| ^_,+ d_d < C||a|| d.^lflfs if s>l-min(|,^) and 1 < p < 4. 

R ^ p T^P P’l ^ 


P,1 
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In order to prove (I102p for /, it suffices to bound {auY in B 2 i)- Now, using 

Bony’s decomposition and the fact that a = + a^, we see that 

(105) {auY = {TauY + {R{a, u)Y + (r„a^)^ + 

The first three terms may be bounded thanks to Prop. 12.31 and Inequalities (11031) . (|104ll 
with s = ^ — 1. Observing that for any Besov norm, we get 


||(i?(a,'u)> 






\^p,l ) 

< C||a|| d_-^ I 


t '■ V 


\u\ 


^ ^““1 


°(S^(oo) 




^+1 ) 
p+) 




LUBl 


^+1)' 


. 4-1 ._ 

Because is embedded in B^^o^ the above right-hand sides may be bounded by 

CXp{t). To handle the last term of (jl05|) . we just have to observe that owing to the 
spectral cut-off, there exists a universal integer Nq so that 


{Tua^^ = Sk,+i ( Sk-iu Aka^) • 

\k—ko\<No 


Hence ||T„a'*|| . d ps 2^“2 J2\k-ko\<No \\Sk-iu ^ka'^Wi^- Now, if 2 < p < m.m{d,2d/{d - 2)) 

Bi,i 

then we may use for \k — ko\ < Nq 

2'^°Y\Sk-iuAkaY\L^ < C2"° WSk-iu^d {2’^i\\AkaY\Lp), 
and if d < p < 4 then 

2'^°Y\Sk-iuAkaY\L^ < C'2'=°(2'=(f-i)||5fc_in||Lp)(2'^p||Afca'*|Up). 

Hence one may conclude that / satisfies (I102p . Bounding g is similar (see |18]1. 


Last step: Global estimate. Putting all the previous estimates together, we get 

(106) Xp{t)<C{Xp{0)+X^{t)). 

Now it is clear that as long as 

(107) 2CXp{t) < 1, 

Inequality pi06p ensures that 

(108) Xp{t) < 2CXp(0). 

Using a bootstrap argument, one may conclude that if -^p(O) is small enough then (f92P and 
(11071) are satisfied as long as the solution exists. Hence (11081) holds globally in time. 
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4.2.1. The proof of Theorem \4.1\ We just give the important steps. We fix some initial 
data so that Xq is small enough. First, Theorem 13.11 implies that there exists a unique 

. ^-1 

maximal solution (a, w) to (f26]i on some time interval [0, r*[, with a E C{f),T*)] 3^^ ), 

l|a||Loo(o,r*xR'i) < 1/2 and u E C{[0,T*); n Ll^^{0,T*; From dUD and Propo¬ 

sition HU one may check that the additional low frequency information is preserved on 
[0,r*): we have 


- 1 , 


eCi[0,T*y,Bi, )nL^(0,r 


-1, 


iii|+') and n'EC([0,r*);il|/)nLL(0,r*;il|r). 


I+F 


Let us assume (by contradiction) that T* < oo. Then applying (llOSp for all t <T* yields 


ll^ll~ — Tllrill^ d.—i — CXq. 

If Xq is so small as (llOSp to imply that both ([3Tp and (|67P are fulfilled on [0, T*) then, 
for all to E [0,r*), one can solve (f26P starting with data {a{to), u(to)) at time t = to and 
get a solution according to Theorem 13.II on the interval [to,T + to] with T independent of 
to. Choosing to > T* — T thus shows that the solution can be continued beyond T*, a 
contradiction. □ 


5. Asymptotic results 

In this section, we focus on two types of asymptotic issues for small global solutions to 
(|Tp that received a lot of attention since the eighties: the low Mach number asymptotic, 
and the long time behavior. We shall see that essentially optimal results may be obtained 
by very simple arguments from the global result we established in the previous section. 


5.1. The low Mach number limit. This subsection is devoted to the rigorous justifica¬ 
tion of the convergence of ([4]) to the incompressible Navier-Stokes equations 


(109) 


dtu + u ■ Xu — fj-Au + Vn = 0, 
divtt = 0, 


in the so-called ill-prepared data case, where we only assume that e ^{po ~ 1) 

suitably bounded. In particular, if we set '=^ — 1); this means that {dta^, dtu^)\t=o 

is of order 1/e, and that one cannot exclude highly oscillating acoustic waves. More con¬ 
cretely, we have to pass to the limit e —)• 0 in: 


( 110 ) 


„ , divu^ , , ,, 

dta -= —div(a M ), 

OF F F Va^ 

dtu + u ■ Xu — — -^ -b 


k{£a^ 


-Xa^ 


1 -b ea^ £ £ 

- ^ - div (2'jl(£a^)D(u^) -b A(ea'^)divu^Id). 

1 -b ea^ 


Before stating our main results, let us introduce some notation. In this section, we agree 
that for z E 5'(M'^), 


( 111 ) 


2i/3<2J0 


^h,l3 d|f 




2J/3>2J0 
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for some large enough nonnegative integer jo depending only on p, d, and on the functions 


def 


k, A/u, /i/u with u = A + 2/x. The corresponding “truncated” semi-norms are defined by 

( 112 ) 


- IP 

p,r 


ana - IP lls^^- 


Keeping in mind the linear analysis we performed for (j70ll in the case u = 1 and e = 1, and 
combining with the change of variable 


(113) 


def 


{a,u){t,x) = e{a^ ,u^){e^ut,ei'x), 


,def 


we expect the threshold between low and high frequencies to be at 1/e with e = eu, and 
it is thus natural to consider families of data (oqjMq) such that 


(oq,mo)|| ’“d + IIoqII 


h,£ I II £1 

d +IpSI 

BP, 


h^e 
. ^-1 
pA 


is bounded independently of e. We expect the corresponding solutions of (|110l) to be uni¬ 
formly in the space defined by 


-T _ . 1 ,^ 


(a^-^ G a(M+; ) n ), 

G a(R+; i?|i) n L^{R+, i/|i), 
n^’^G a(M+;i?|”') n Li(M+;i?|^\ 


and endowed with the norm: 


[a,u)\\xi„ = 


def,,, 

“>^)IP d 




+ 11^ 


I 

+u\\{a,u)\\ 


+ e a 


I 


£,e , II \\h,e , —In \\h,e 

d,i d,i +£ l|o|| d ■ 

LHBl, ) lHBL ) LHBl,) 


One can now state our main result of convergence in the small data case, the reader being 
referred to da El for the large data case and stronger results of convergence. 


Theorem 5.1. Assume that the fluid domain is either Mfl or that the initial data (oq, Mq) 
are as above and that p is as in Theorem There exists a constant p independent of £ 
and of V such that if 


(114) 


^£,!^ I 


(«0>«0; 


^2,1 


+ 


+^l“oll^’d <w, 

BP, BP, 

p, 1 p, 1 


then System (IllOjl with initial data (ogjUg) has a unique global solution {a^,u^) in the space 
with, for some constant C independent of e and v, 


(115) 


lx?,. < CCIX 


In addition, Qu^ converges weakly to 0 when e goes to 0, and, if Vuo ^ vq then Vu^ 
converges in the sense of distributions to the unique solution of (jlOQh supplemented with 
initial data vq- 


Proof. Performing the change of unknowns given in ()113p and the change of data 
(116) {ao,uo){x)^= £{al,ul){£i^x) 
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reduces the proof of the global existence to the case u = 1 and e = 1, which was done in 
Theorem KT\ Back to the original variables will yield the desired uniform estimate (11151) 
under Condition (jll4p . Indeed, we notice that we have up to some harmless constant: 


4^1 


+ 


e\\h,e 


. ^-1 
BP, 


+ £ a) 


oil d 

BP, 

p,i 


= z^(||(ao,^io)|| 




^2,1 


I II ll^d I II I 

+ IPoII + ||flo| 


BP, 

pA 


h,l 

d 

pA 


) 


and 




Granted with the uniform estimates established in the previous section, it is now easy to 
pass to the limit in the system in the sense of distributions, by adapting the compactness 
arguments of P.-L. Lions and N. Masmoudi in |24j . 

More precisely, consider a family (aQ,nQ) of data satisfying (|114ll and Vuq vq when 
e goes to 0. Let {a‘^,u^) be the corresponding solution of (IllOl) given by the first part of 
Theorem O Because 


(117) 




BP, 

pA 






the data (oq, Uq) are uniformly bounded in ^ i- Likewise, (11151) ensures that {a^,u^) 

■ 4-1 

is bounded in the space Cb(]R+; ). Therefore there exists a sequence (en)neN decaying 

. 4-1 

to 0 so that ^ (ao,Mo) in ^pi (with Vuq = uq) and 

. 4-1 

(118) ,u^") ^ {a,u) in L°°(R+-,Bp^ ) weak *. 

The strong convergence of the density to 1 is obvious: we have = l+enO^", and (0'^")nGN 

. 4 

is bounded in L^(M+; Bpi)- In order to justify that divu = 0, we rewrite the mass equation 
as follows: 

divu^" = — endiv(a^"u^") — Sndta^". 

d 

Given that and are bounded in L^(M+; Bp i) (use the definition of and interpo- 

• 4—1 

lation), the first term in the right-hand side is 0{sn) in L^(]R+; Bp^ ). As for the last term. 


it tends to 0 in the sense of distributions, by virtue of (|1I8I) . We thus have divu^ 
whence divu = 0. 


0 , 


To establish that u is a solution to (II09p . let us project the velocity equation onto 
divergence-free vector fields: 


(119) 


dtVu^" - = -P(n^" • Vu^") - V 


1 


-Au^ 


, 1 -7 

Because Qu = 0, the left-hand side weakly converges to dtu — fiAu. To prove that the last 

d 

term tends to 0, we use the fact that having and (a^)^’^ bounded in L°°{Bp^) and 


■ 


L°°{BP,^ ), respectively, implies that, for all a E [0,1], 




■ --l+a 

is bounded in L°°{Bp -^ ). 


(120) 
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■ --1 

Now, Au^ is bounded in ) and p < 2d. Hence, according to product laws in Besov 

^ . ^—2+a 

spaces, composition inequality and ()120p . we get {l+£a‘^)~^Au^ = in ), 

whenever 2max(0,1 — ^) < a < 1. Hence the last term of pi9p goes strongly to 0 for some 
appropriate norm. 

In order to prove that ^ T’{u ■ Vu), we note that 

Projecting the first term onto divergence free vector helds gives 0, and we also know that 
Vu = u. Hence we just have to prove that 

(121) 'P{Vu^'^ ■ V{Vu ■ Vu) and V{Qu^'^ ■ ^ 0. 

This requires our proving results of strong convergence for Vu^’^. To this end, we shall 
exhibit uniform bounds for dt'Pu^'^ in a suitable space. First, arguing by interpolation, 

^ _|_ 2 _g 

we see that is bounded in ) for any m > 1. Choosing m > 1 so 

that ^ — 3 > —(imin(|, 1) (this is possible as p < 2d) and remembering that (e^a^") is 

bounded in L°°{B^^), we thus get {{l+£na^"")~^Au‘^'^) bounded in ). Similarly, 

. ^-1 . 

combining the facts that (tt^") and (Vn^") are bounded in L°°(Hp ) and L'^{Bp^ "" ), 

d _|_ 2 _ 2 

respectively, we see that ■ Vu^") is bounded in L'^{Bp^ ™ ), too. Computing dtVu^" 

from (|119p . it is now clear that {dtVu^") is bounded in L™'(i?p ^ ). Hence {'Pu^'^ —Vu^") 

is bounded in (M_,_; ). As Vu^’^ is also bounded in Cb{^+]B^i ), and as the 

embedding of ^ in B^^ is locally compact (see e.g. [T], page 108), we conclude by 
means of Ascoli theorem that, up to a new extraction, for all (f) G 5(M'^) and T > 0, 

(122) ^Vu^'^^4>Vu in C([0 ,T];h|J'™"% 

Interpolating with the bounds in C;,(M+; B^^ ), we can upgrade the strong convergence in 

■ - —1 —Q 

(11221) to the space C([0, T]-,Bp^ ) for all small enough a > 0, and all T > 0. Combining 

with the properties of weak convergence for Vu*"" to Vu, and to 0 that may be 

deduced from the bounds of it is now easy to conclude to (jl2ip . One can use for 
instance the fact that for all m > 1, we have 

^ Vu in L^(B^^^~^) weak * and ^ 0 in L™(H|^™"^) weak * . 

□ 

5.2. Time decay estimates. In the present subsection, we show that under a mild addi¬ 
tional decay assumption that is satisfied if the data are in L^(R'^) for instance, the norm 
(the i ?2 1 norm in fact) of the global solutions constructed in Theorem 14.11 decays like t~i 
for t —)• -boo, exactly as for the linearized equations. This fact has been hrst observed by 
A. Matsumura and T. Nishida in [25j in the case of solutions with high Sobolev regularity. 
The adaptation to the type critical regularity framework has been carried out recently 
by M. Okita in [27], in dimension d >3. Below, we give a more accurate description of the 
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time decay, emphasizing a better decay for high frequencies. This is the key to handling 
any dimension d >2. For simplicity, we concentrate on the type framework, even though 
we expect similar results to be true in the more general framework of Theorem 14.11 

Theorem 5.2. Let the data (ao,uo) satisfy the assumptions of Theorem \f.l\ with p = 2 and 

assume with no loss of generality that P'{1) = 1 and that v = 1. Denote (r) ='^ vT+T^ 

and a = min(^ + 2, | + ^ — e) with e > 0 arbitrarily small. There exists a positive constant 
c so that if in addition 

(123) Do sup (||.T(Afcao)||Loc + ||.F(AfcUo)||Lco) < c 

k<ko 

then the global solution {a,u) given by Theorem \4.1\ satisfies for all t>0, 

(124) D{t) < C{Do + ||(oo, Vao,uo)|| . d.J 


def 


with D{t) = sup ||(r)4 + 2(a,u)|| + ||(T)“(Va,u)|L . <i_i + ||rVu|| 




sG(-f,2] 

Proof. Throughout the proof, we shall use repeatedly that for 0 < cii < cj 2 , we have: 

(125) f {t — t)~^^ dr < {t)~^^ if in addition (T 2 > 1. 

Jo 

Step 1: Bounds for the low frequencies. Denoting by E{D) the semi-group associated to 
dZQD, we have for all A: G Z, 


(126) 


Afca(t) \ ^ ^tE{D) f AfcOo ^ _ 

Afcu(t) J V j^kUO 


,(t-T)E{D) 


Afc/i(r) 

Afc(/2 + /3 + /4)(r) 


dr 


with fi div(an), f 2 ’=^ u ■ Vn, /a k{a)Va and fi /(a).4.u. 

From an explicit computation of the action of in Fourier variables (see e.g. m, 

we discover that there exist positive constants c and C depending only on ko and such that 

\P{e^^^^^U){f)\ < for all |C| < 

Therefore, for all k < ko, 

< ||.FAfcC/||i.o2“e-'=°22''*. 

We thus get up to a change of cq, 

(127) t^+f ^ 2'=^||e*^(^)Afc[/|A2 < (sup WPAkUU^) ^ (Vt 2 *=)^+'* 

k<ko k<ko 

As for any cj > 0 there exists a constant so that 

(128) sup 

fcez 

we get from (I127p that for s > —d/2, 

sup < a sup ||.FAfct/|Aoc. 

i >0 2.1 
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It is also obvious that for s > —df2, 

||gtE(D)^||r 


L < 


Bs < sup ll-FAfeC/llio 
'“2.1 k<ko 


Hence we conclude that 


(129) 


sup(t)i+§||e'®(^)[/||<._ < sup \\J^AkU\\B 


t >0 2,1 

Next, we claim that for all s E {—d/2,2] and i E {1, • • • ,4}, we have 


(130) 


/ (t-r)" 3 -f sup \\TAkfi{T)\\L^ dr < + X^{t)) 

Jo k<kQ 


with X(t) ||(a, Va,u)||_ d_^ + [ (||a||^ d , ^ + ||a||^ d + llitll . 4 + 1 ) dr- 

') 7o ' Bit Bl, Blt^' 

Of course, as the Fourier transform maps to L°°, it suffices to prove (jl30p with ||/j|| 2 ,i 

instead of sup^^^^, \\TAkfi\\L^. 

To bound the term with /i, we use the following decomposition: 

fi = u ■ Va + adivu^ + adivu^. 

Now, from Cauchy-Schwarz inequality and the definition of D{t), one may write 

[ - r)“^“^||(u • Va)(r)||ii dr < ( sup (r)^||u(r)||i 2 ') ( sup (T)5+i||Va(r)||i 2 ') 

Jo ^0<T<t 2 \o<r<t 2 

X dT 

Jo 


where we used (jl25p and the fact that 0<^ + |<| + |- 

Bounding the term with adivu^ is totally similar. Regarding the term with adivu^, we 
use that if t > 2, 

[ {t-r)~i~^\{adivu^){T)\\LidT < j (t - T)“3-f ||a(r)||i 2 ||divu'"(r)||i 2 dr 
Jo Jo 

Therefore, as —d/2 < s < 2, we get 

(t)t+3 f (t - r)"^“5||(adivu'*)(r)||£i dr < ( sup ||a(r)||i 2 ') f ||divM''(r)||i 2 dr 
do Ve[o,t] 2 Jq 

+ ( sup (r)^||a(r)||i 2 ) ( sup r||divM^(r)||i 2 ), 

VG[0,t] / 


Tg[ 0 ,t] 


and (jl30p is thus satisfied by the term with adivu^ if t > 2, the case t < 2 being obvious 
as (t) ~ 1 and (t — r) ~ 1 for 0 < t < t < 2 and one may write 

adivn^ll^i dr < ||a|L 2 (^ 2 )||divu'^||i 2 (^ 2 ) < X'^{t). 
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Handling the terms with /2 and /a is totally similar: k{a)S/a and u ■ may be treated 
as u ■ Va, and u ■ Vu^, as adivu^. For f^, we write that 

/4 = I{a)Au^ + I{a)Au^. 


Now, we have 

rt 


f {t — t) 4 2 \\I(^a)Au^\\Li dr 

Jo 

< ( sup {T)i\\a{T)\\L 2 ')( sup (r)^+^||V^M^(r)||i 2 ') f {t - dr. 

VG[0,t] ^ ^ Jo 

Hence, thanks to (I125h . the term with I{a)Au^ fulfills (I130p . Finally, for t >2, 

[ (i -dr < / ||a||2,2 1 |dr 

do do 

+ ^ ((r)i ||a(r)||i2) (r||V2M^(r)||i2)dr, 


hence, because —d/2 < s <2 and < ||rVrt|| 




f 


(t-r) 4 2 ||/(a)^tt^||^i dr < (t) 4 2 (^d"^( t) + X‘^{t)) for t>2. 
Obviously, as (t) 1 and (t — r) « 1 for 0 < r < t < 2, we have the following inequality: 


f 


{t — t) 4 2 \\I(^a)Au^\\Li dr < (t) 4 2X^{t) for t<2, 


which completes the proof of (|130|) . Combining with (I129p and using Duhamel’s formula, 
we conclude that for all t > 0 and s G (—d/ 2 , 2 ], 


(131) 


(t)4 + 2||(a,u) 


\bs <Do + D\t)+X\t). 

JDO T 


Step 2: Decay estimates for the high frequencies of (Va,u). We here want to bound the 
second term of D{t). Recall that Theorem 14.11 ensures that 

||(Va, 7 r)|l <CX(0) for all T > 0. 

^t(^2,1 ) 

Therefore it suffices to bound |h"(Va,ri)|l d_^ for, say, T >2. 

(^ 2,1 ) 

Now, the starting point is Inequality (j87p which implies that for k > and for some 
Co = c{ko) > 0 , we have 

^±Cl + coCl<{\\{Vh,g,)h2 

+ \\Rk{u,a)\\]^2 + \\Rkiu,u)\\i2 + \\Rk{u,a)\\i2 + \\Vu\\loo£i.'^£i. 
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After time integration, we discover that 

<£fc(0) + £e^“"(||(V/fc,fffc)||i2 + ||i?fc(n,a)||i2 

+ \\Rk{u,u)\\L2 + \\Rk{u,a)\\L2 + \\Vu\\L^Ck^dT, 

whence, remembering that Ck ~ \\{AkVa, Aku)\\i 2 for k > /cq, 

II (AfcVa, Aku) (t) 11^2 < rII (AfcVa, Afcn)(0) 11^2 + (|| (V/fc, 5fc) Hl^ 

+ ||i?A:(^^,a)||L2 + \\Rkiu,u)\\L2 + ||^fc(u,a)||i2 + llVullLoclKAfeVa, AfcM)||i2^<ir 

and thus, multiplying both sides by taking the supremum on [0,T], and summing 

up over k > ko, 

(132) ||r(Va,n)L . < ll(Vao, uo)f sup dr) 


with Sk 

Sl=^mfk,gk)\\L^, Sl^^^\\Rk{u,a)\\L2, \\Rk{u,u)h2, 

St‘^^^\\Rk{n,a)h2, tf||Vu|U^||(AfcVa,Afen)|A2. 

In order to bound the sum, we first notice that 



^o(^-*)2^(5-i)5fc(r)dr^ 



k>ko 


Hence taking advantage of ()24ll and of a similar inequality for Rk{u,a), we end up with 

^ sup t“/W-*)2"(i-i)5fcdr< r(||(V/, 5 )||.|_i + ||Vu||.. ||(a, Va, n)|| . ._,)dr. 
I^^0<t<2 Jo Jo ^ Bl^ Bl^ / 


Bounding V/ and g as in the proof of Theorem 14.11 leads to 

(133) V sup t" / < a2(2). 

rri o<t<2 Jo 


k>kQ — — 


To bound the supremum on [2,T], we split the integral on [0,t] into integrals on [0,1] and 
[l,t], respectively. The [0,1] part of the integral is easy to handle: we have 

V sup r Skir) dr < V sup 2^(i-^)5fc dr 

^k,^<t<T Jo fctt2<t<T Jo 

k>kQ 


Hence 

( 134 ) 


E 

k>ko 


sup 

2<t<T 
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Let US finally consider the [l,t] part of the integral for 2 <t <T. We shall use repeatedly 
the following inequality 


(135) 


llrVnll 




< 


D{t), 


which is straightforward as regards the high frequencies of u and stems from 


||rVu|| 


1 4 <||(T)^Wu||i . <|Kr)l+^u| 




<m 


for the low frequencies of u. 

Regarding the contribution of 5^, we first notice that, by virtue of (I125j) . 

1 I e -Okyijwi - ||7 y\/j,yj\'^ 

Now, product laws in tilde spaces ensures that 


(136) V sup < ||r“(V/,, 

7^ 2<t<T Jl 


4-1 • 


k“v/|| 


- 4-1 ^ lk“ d ||rdivM||„ d . 


The high frequencies of the first term of the r.h.s. is obviously bounded by D{T). As for 
the low frequencies, we notice that if d < 4 then for all small enough e > 0, 


I ^-e I 

r2 a 


^ II 4 _ £ 

< lira ^a l 


(137) „ _ 

and if d > 5, taking s = 2 in the first term of D{T), 

(138) 




< D{T) 




Therefore, using ()135l) and remembering the definition of a, we get 




Next, we have 


|r"(A:(a)Va'^)|| 4_, < ||a| 


11^ 

d T CL\\ 


_ , < X{T)D{T) 


and, according to (I137|) and 


|T“(A:(a)Va^)L ^.4 <||ri--a| 

BriB^y ) 


d ||r"-i+^a|K 


< 


L5?(R^i) 




D‘^{T). 


We also see that 
||r"/(a)^u| 




-1 ^ 


|r"-^a|K 


B¥(B2,1 ) 


I ||T““^nll^ 

_ d W ||7 ^IL 4 
B¥(B2,i) 


)• 


The first term of the r.h.s. may be bounded by virtue of (11351) . and it is also clear that 
the last term is bounded by D{T). As for the second one, we use again ()137p and (|138p . 
Resuming to (jl36p . we end up with 

V sup r/ e'^°(^-*)2^(i-i),S^(r)dr < Z12(T). 

fc^„2<t<T Jl 
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To bound the term with we use the fact that 


^'^\\Rk{u,a)\\L 2 dT <\\Rk{TU,T°‘ ^a)||ioo(^ 2 ) ^ “dr. 

Hence, thanks to ()125l) and to (l2^ (adapted to tilde spaces), 
sup (t°‘ f ^ 

^.■ 2 <t<T\ Ji / 


k^ko — — 


k^ko 

||rVu|| 


d ||r“ ^a|L d, . 


The hrst term of the r.h.s. may be bounded thanks to (11351) . and the high frequencies 
of the last one are obviously bounded by D{T). To bound ||r““^o||^ , we use the 

following two inequalities 


|r"-^a|K 


~ • w — 1 

l5?(B2 i ) 


T a 




I «—1 id 

r a 






lL5?(B| i) 


if d < 6 , 
if d > 7. 


^ _ 1 
2 2 

by D{T). We eventually get 


Because a — 1=2 — 5 — eifd<6, and a — l = | + lifd>7, the r.h.s. above are bounded 


Y sup r < Z12(T). 




The terms and may be treated along the same lines. 
Finally, using product laws and (|125l) . we get 


V sup 

2<t<T Jl 


k>ko — — 


< 


||rVu||„ ||r“-i(Va,u)||(( sup dr < D‘^{T) 


Lt(^2,i) 


) 2<t<'r Jl 


Putting all the above inequalities together, we conclude that 

V sup (r [ <D{T)X{T) + D^{T). 

k>k,^<t<T\ Jl ) 

Then plugging this latter inequality, (I133h and (I134h in (I132h yields 

_ ^.d. < ||(Vao,uo)||%_, +x2(r) + di2(r). 


(139) 


(T)“(Va,u)|| 


Step 3: Decay estimates with gain of regularity for the high frequencies of Xu. In order to 
bound the last term of D{t), we shall use the fact that the velocity u satisfies 

Hpf 

dtu — Au = F = —(1 + k{a))Xa — u ■ Xu — I{a)Au, 

whence 

dt{tAu) — A{tAu) = Au + tAF. 
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Because the maximal regularity estimates for the Lame semi-group are the same as for the 
heat semi-group (see the beginning of Section [3|), we deduce from Remark 12.21 that 




d 1 ^ ll'4.u||^ d_^ + .d_3 


^_1 

BriBL ) 


whence, using the bounds given by Theorem 15.1 


(140) 


||tVtt| 


LTiBh) 


<X(0) + ||rFf 


LTi^A S 


In order to bound the last term, we notice that, because a > 1, we have 


||rVa|| 


" .d , <||(r)“a||'( ,. 


Next, product and composition estimates adapted to tilde spaces give 


as well as 


and 


|r A:(a)Va||'( d < Il'r^aH^ d < 


LTiBli 


|t n • Vn||( 


< 


) 


1 d_i llrVulL d 

LriBl^ ) LT{Bh) 


|r/(a).Au||_ d 


d , < 

LriBi^ ) 

Therefore, resuming to (11401) and remembering (11351) . we get 


1 ~ II® ~ d rV u L d_, 

^ ' LTiBh) ) 


\0'^\\h 


. <X( 0 ) + Zl(t)X(t) + Zl^(t) + ||(r)“a|L d 

LriBli) L^iBl,) 


Finally, bounding the last term according to (I139p . and adding up the obtained inequality 
to (I13ip and ()139p yields 


Dit) <Do + IKVao,no)f , , + X^t) + D^it). 


B^~ 

^2,1 


As Theorem 15.11 ensures that X{t) is small, on can now conclude that ()124l) is fulfilled for 
all time if Dq and ||(Vao,no)||^ d_, are small enough. □ 


^2,1 


6. Appendix 

Here we recall various estimates for the flow that have been used repeatedly in the proof 
of Theorem EH More details may be found in m or m- 

Recall that if v : [0,T) x —)• is measurable, such that 1 1 -)> v(t,x) is in L^(0,T) for 
all X € R*^ and in addition Vu G L^(0, T; L°°) then it has, by virtue of the Cauchy-Lipschitz 
theorem, a unique flow X^ satisfying 

Xv{t,y)=y+ [ v{T,X^{T,y))dT for all tG[0,r). 

Jo 

In addition, for all t G [0,T), the map Xy{t, •) is a C^-diffeomorphism over R'^. 
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(IqJ' 

Lemma 6.1. Let p ^ [IjOo). Letv{t,y) = v{t, X{t,y)). Under Assumption (l60l) . we have 
for all t G [0, T], 


(141) 

(142) 

(143) 

(144) 


\\ld-adi{DX,m\,a <||I11;|| 






I|ld-^4i)|| . <||Z11;|| 


<1 


LliBU) 


||adj(IlX,(t))'^^,(t)-Id|| < \\Dv\\ 


RP 


LliBfi) 




< 


BP, 

p.i 


\\Dv\\ 


LliBfU 

Proof. As an example, let us prove the last item. We have thanks to the chain rule, 

(145) Jjj{t,y) = l+f div v{T,Xjj{T,y))Jjj{T,y) dr = 1+ f (Du : adj (DX^))(r, y) dr. 

Jo Jo 

Hence, if Condition (1601) holds then we have (I144p for J„, a consequence of the fact that 


Bp ^ is an algebra, and of (I14ip . In order to get the inequality for ^, it suffices to notice 
that 

- 1 = (1 + {Jv{t,y) - 1))"^ - 1 = [ Dv : adj {DX^)dT. 


k>l 


□ 


def _ 


Lemma 6.2. Let vi and V 2 be two vector-fields satisfying (l60p . and 6v = V 2 — vi. Then we 
have for all p G [1, oo) and t G [0, T] .• 


(146) 


14 — 4 II 

\-^V2 ^Plll .= ^ 

LTiBU) 


<||D(&;|| , 

LiiBfU 


(147) 


11 adj {DXy^) - adj {DXy 


< 


LpiBP) 


\\D6v\\ 


. d ) 
Ll(BU) 


lHbL) 


(148) 11^.2-^.1II 4 <ll^<5u|| 

Proof. In order to prove the first inequality, we use the fact that, for i = 1,2, we have 

Ay, = (Id + Ci)-^ = V(-l)^Cf with Ciit) = f Dvi dr. 

k>o -^0 


Hence 

H,, - H,, = - Cf) = r fodvdA 

k>l ^ fc>lj=0 

. ^ 

So using the fact that ^ is a Banach algebra, it is easy to conclude to (114611 . Proving 
the second inequality is similar. □ 
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